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(57) Abstract: The invention relates 
to a method for the electromagnetic 
modification of cerebral activity in order 
to reach a pre-determined behaviour target, 
using a behaviour model which describes 
the correlation between behaviour and the 
course of cerebral activity characteristics, 
JQ and a cerebral activity model which 

quantitatively describes cerebral activity. 
The cerebral activity characteristics and 
the courses thereof can be derived from 
the quantitative description of the cerebral 
activity. According to the inventive method, 
the nominal course of given cerebral activity 
characteristics which corresponds to the 
achievement of the behaviour target is 
determined by means of the behaviour 
model, exogenous electrical and/or magnetic 
fields are produced, and electromagnetic 
cerebral activity is measured and cerebral 
activity characteristics are calculated in 
order to achieve the nominal course in a 
control or regulating circuit containing 
calculations with the cerebral model. 
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(57) Zusammenfassung: Die Erfindung betrifft ein Verfahren zur elektromagnetischen Modifikation von Himakti vital zum Errei- 
chen eines vorgegebenen Verhaltenszieles, unter Verwendung eines den Zusammenhang zwischen Verhalten und den Verlauf von 
Himaktivitatscharakteristika beschreibenden Verhaltensmodells, sowie eines Hirnaktivitatsmodells, welches Hirnaktivitat quantita- 
tiv beschreibt, wobei Himaktivitatscharakteristika und deren Verlaufe aus der quantitativen Beschreibung der Hirnaktivitat ableitbar 
sind, umfassend Ermitteln mittels des Verhaltensmodells, welcher Soll-Verlauf bestimmter Hirnaktivitatscharakteristika dem Er- 
reichen des Verhaltenszieles entspricht, Erzeugung exogener elektrischer und/oder magnetischer Felder, und Messung elektromag- 
netischer Hirnaktivitat nebst Berechnung von Hirnaktivitatscharakteristika zum Erreichen des Soll-Verlaufs in einem Steuer- oder 
Regelkreis, der Berechnungen mit dem Hirnmodell beinhaltet. 
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VERFAHREN UND VORRICHTUNG ZUR ELEKTROMAGNET1SCHEN 
MODIFIKATION VON HIRNAKTIVITAT 

Die Erfindung betrifft ein Verfahren und eine Vorrichtung zur elektromagnetischen 
5 Modifikation von Hirnaktivitat, insbesondere fur die gesteuerte oder geregelte 
elektromagnetische Modifikation von Hirnaktivitat in vivo und die daraus 
resultierende Verhaltensmodifikation. 

Definitionen: 

10 

Modellbasierte gesteuerte oder geregelte Modifikation bedeutet gezielte 
Veranderung oder Nichtveranderung von Hirnaktivitat auf der Basis eines 
Verhaltenszieles, eines Verhaltensmodells und eines Himaktivitatsmodelis. 

15 Geregelt bedeutet, dass Beobachtungs-, Berechnungs- und Modifikationsschritte 
in Ruckkopplungsschleifen miteinander verbunden sind. 

Das vom Anwender zu benennende Verhaltensziel umfasst Art und Auspragung 
eines bestimmten Verhaltens sowie dessen Dauerhaftigkeit. 

20 

Verhalten bezieht sich auf jegliche spatiotemporale Kombination jeglicher 
potentiell validierbarer faktischer oder moglicher Zustande einer Person. 
Verhalten beinhaltet unter anderem, aber nicht ausschlielilich, Wahmehmung, 
Handlungen, Motivation, Aufmerksamkeit, Gedachtnis, Lemen, Bewusstsein, 
25 Emotionen, Denken, Kognition, mentale Zustande, Fahigkeiten, 
Personlichkeitseigenschaften. 

Verhalten basiert auf Hirnaktivitat. Der Zusammenhang zwischen beiden wird im 
Rahmen eines Verhaltensmodells quantifiziert, d.h. einer empirisch festgestellten 
30 Verknupfung von bestimmtem Verhalten mit einer bestimmten Dynamik 
bestimmter Hirnaktivitatscharakteristika. 

Hirnaktivitatscharakteristika sind aus Messwerten ermittelbare Kenngrofien von 
Hirnaktivitat. Beispiele fur Hirnaktivitatscharakteristika sind Potentialdifferenzen 
35 zwischen einer EEG-Elektrode und einer Referenzelektrode, oder der 
Bandbereich der uber alle EEG-Elektroden ausgepragtesten Frequenzen („Alpha- 
Band" oder ahnliches [12]), oder der Ahnlichkeitsindex (^similarity index", siehe 
[5]) der Messwerte von MEG-Sensoren. 



40 Ein Hirnaktivitatsmodell bezieht sich auf ein physiologisch basiertes Modeil, 
welches Hirnaktivitat vermittels der Dynamik von Hirnbestandteilen und/oder 
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deren Wechselwirkungen beschreibt (das umfasst z.B. in [24], [1], [23], [25] oder 
[18]). 

Konventionelle Modelle sind nicht spezifisch, sondern generisch , da 
5 Nichtobservable uber Populationsmittelwerte geschatzt werden. 

Eine Nichtobservable ist ein Parameter Oder eine Variable eines generischen 
Modells, der/die in vivo nur mit erheblichem Aiifwand - oder gar nicht - mit der 
notigen raumlichen und zeitlichen Auflosung gemessen werden kann. 
10 Nichtobservable bestimmen entscheidend die Hirndynamik, und schwanken 
zwischen Individuen und innerhalb eines Individuums bis zu einem Faktor von 
Einhundert. Generische Himaktivitatsmodelle machen keine Aussagen uber 
direkte Fremdeinwirkung auf die Dynamik von Hirnbestandteilen, bestenfalls uber 
indirekte Fremdeinwirkung durch Rezeptorstimulation, Deprivation oder Lasionen. 

15 

Die im folgenden verwendeten, im Regelfall spezifischen Himaktivitatsmodelle 
quantifizieren unter anderem Fremdeinwirkung uber elektrische und/oder 
magnetische Felder und basieren auf fur einen Anwender in einem Zeitintervall 
relevanten, mit Hilfe geeigneter Verfahren ermittelten Nichtobservablen. Diese 
20 Himaktivitatsmodelle umfassen fur jeden Himbestandteil in der Regel nichtlineare 
Differentialgleichungen beispielsweise der Gestalt 

x(t) = f(xjJokaleParameter,endogenerInput,exogenerInpu^ 

25 Hierbei ist x ein Himaktivitatscharakteristikum, x dessen Zeitableitung, t die Zeit, 
f eine Funktion unter anderem der Art des Bestandteils (Jokale Parameter"), der 
Art und Umfang der sich uber physiologische Verbindungen ausbreitenden 
Einwirkung anderer Bestandteile („endogener Input 41 ), und der direkten 
Fremdeinwirkung auf den Himbestandteil („exogener Input"). Der endogene Input 

30 umfasst abhangig vom verwendeten Modell weitere („translokale") Parameter, 
wie bspw. Kopplungsstarken oder Verzogerungszeiten zwischen Bestandteilen. 
Fehlende oder unzureichende Spezifikatlon des Hirnaktivitatsmodells auf den 
jeweiligen Anwender fuhrt zu einer breiten Streuung in diesem Fall weitgehend 
unbekannter Auswirkungen exogenen elektromagnetischen Inputs und somit zur 

35 Nichtanwendbarkeit entsprechender Verfahren aulierhalb eines Bereichs 
neurologischer, psychiatrischer, oder psychologischer Schadigungen, in dem 
Zufallsauswirkungen gegen mogiichen Behandlungserfolg abgewogen werden. 
Elektromagnetisch bedeutet elektrische und/oder magnetische GrofJen und 
bezieht sich sowohl auf die Art der eingehenden Beobachtungsdaten (im 

40 wesentlichen elektromagnetische Korrelate der Hirnaktivitat des Anwenders), als 
auch auf den wesentlichen Modus der Modifikation (beispielsweise per 
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extrakranial erzeugtem veranderiichen Magnetfeld, welches miihelos das 
Kranium des Anwenders durchdringt und intracranial Induktionsspannungen 
hervorruft, wodurch eine Einflussmoglichkeit auf die Hirnaktivitat des Anwenders 
gegeben ist). 

5 

In vivo bezieht sich auf die Verwendung des Verfahrens bei lebenden 
Anwendern, im Gegensatz zu Zellpraparaten, Hirnteilen, Computersimulationen 
oder ahnlichem, und stellt somit besondere Anforderungen sowohl im Hinblick auf 
die Komplexitat der stattfindenden Hirnprozesse, als auch an die Beobachtungs-, 
10 Berechnungs- und Modifikationsgeschwindigkeit sowie Sicherheit des 
Verfahrens. 

Relevante Technologien: 

15 Mit dem gegenwartigen Stand der Technik ist modellbasierte gesteuerte oder 
geregelte elektromagnetische Modifikation von Hirnaktivitat in vivo nicht moglich 
(das menschliche Gehirn besteht aus circa 10 11 Neuronen und circa 10 15 
Verbindungen zwischen diesen), und wurde bis dato weder in Betracht gezogen, 
noch versucht. Das gleiche gilt fur die damit zusammenhangende 

20 Verhaltensmodifikation. 

Es gibt drei, wenige Uberschneidungen aufweisende wissenschaftlich / 
technische Entwicklungsrichtungen mit Bezug zu Hirnaktivitat bzw. deren 
Veranderung: 

25 

einen nichtlineare Dynamik / neurale Netze / kiinstliche Intelligenz - 
Ansatz (im folgenden als "NA" bezeichnet), bei dem bestimmte 
Eigenschaflen des lebenden menschlichen Gehirns in vereinfachter 
Weise nachgebildet werden sollen, urn Maschinen mit ahnlichen 
30 Eigenschaften zu bauen (sog. „Neurocomputer"), 

einen kontrolltheoretischen Ansatz (im folgenden als "CON" 
bezeichnet), in dem Prinzipien der mathematischen Kontrolltheorie auf 
engumschriebene Teile einer Gehirn-Korper-Einheit angewendet 
werden, 

35 - einen medizinischen Ansatz (im folgenden als "MED" bezeichnet), der 
sich mit praktisch-klinischen Verwendungszwecken von 
elektromagnetischen Feldern (im wesentlichen) hoher Feldstarken bei 
Menschen oder Tieren befasst. 

40 Ein moglicher Grundbaustein von NA ist der neurale Oszillator, das heilit eine 
Entitat, die zwischen oszillatorischem Verhalten und Rune wechseln kann, und 
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die sich durchaus aus einem Ensemble kleinerer Entitaten zusamrnensetzen 
kann (die aber nicht weiter betrachtet werden). NA befasst sich mit neuralen 
Netzen, das heifct kleinen (GrolJenordnung < 10 4 ) verbundenen Ensembles 
neuraler Oszillatoren, und den Phanomenen, die in diesen Netzen auftreten (z.B. 
5 Gedachtnis oder Mustererkennung). Nachdem ursprunglich Funktionsprinzipien 
aus Nervensystemen von Lebewesen extrahiert wurden, bewegt sich NA 
zunehmend in Richtung kunstlicher Systeme auf der Basis der extrahierten 
Prinzipien. Exemplarisch fur diese Entwicklung ist [2]. 

10 CON bewegt sich in umgekehrter Richtung, nachdem dieser Ansatz mit der 
Anwendung linearer Kontrolle auf hinreichend einfache mechanische Systeme 
seinen Anfang nahm. Lineare Kontrolle auf nichtlineare respektive stochastische 
Systeme (wie z.B. ein menschliches Gehirn) anzuwenden, ist nicht nur 
ungeeignet, sondern auch potentiell gefahrlich, da bei linearer Kontrolle die 

15 Kontrollkraft proportional zur angestrebten Veranderung ist und somit sehr groB 
werden kann. 

Kontrollverfahren lassen sich in modellbasierte und datenbasierte Verfahren 
einteilen. Die Qualitat eines modellbasierten Verfahrens ist verknupft mit der 
Eignung des Modells fur das in Frage stehende Problem. Bis dato werden im 
20 Humanbereich oder bei Tieren keinerlei Hirnaktivitats- oder Verhaltensmodelle 
eingesetzt. 

Die Qualitat eines datenbasierten Verfahrens ist von der Einfachheit des zu 
kontrollierenden Systems abhangig (da der erste Schritt in einem datenbasierten 
Verfahren in der Regel die Rekonstruktion des Phasenraumes ist, in dem 

25 anschlieliend Suchverfahren erfolgen und Iterationen gerechnet werden - was 
bei hoherdimensionalen Phasenraumen die konkrete Berechnung in Realzeit 
unmoglich macht. Ausnahme: zeitverzogerte Ruckkopplung, bei der jedoch 
Wartezeiten bis zum Erreichen eines zu stabilisierenden Zielzustandes auftreten). 
Fur eine gesteuerte oder geregelte elektromagnetische Modifikation von 

30 Hirnaktivitat in vivo sind gegenwartige Kontrollverfahren nicht anwendbar, da 
jeweils mehrere der folgenden Effekte zusammen kommen: 

hohe bis unerfullbar hohe Anforderungen an Speicherplatz und 
Rechnergeschwindigkeit, 

potentiell beliebig lange Wartezeiten bis zum Erreichen eines Zielorbits, 
35 - intermittierende Ausbriiche aus dem Zielorbit in Systemen mit 
stochastischen Elementen, 

Einschrankungen des Kontrollbereichs in Systemen, in denen nur 
aggregiertes Verhalten messbar ist, 

Einschrankungen des Kontrollbereichs in Systemen, in denen die 
40 Kontrollkraft verschmiert wird, 
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der Betrag der Kontrollkraft im Inneren heterogener Systeme ist 
unbekannt. 

MED beruht implizit auf einem mikroskopischen physiologischen Modell mit 

5 einzelnen Neuronen als Konstituenten der Hirnaktivitat- die im folgenden 
"Einzelneuronen" genannt werden sollen. Demgemafc spielen in vivo nichtinvasiv 
nicht beobachtbare GrolJen wie lonenkonzentrationen auf beiden Seiten der 
Zellmembran, Anzahl der Dentriten etc. eine wesentliche Rolle fur die 
Berechnung der Anzahl der Nervenimpulse pro Zeiteinheit (fire rate). Fur eine 
10 automatische, individualisierte, gesteuerte oder geregelte elektromagnetische 
Modifikation von Hirnaktivitat in vivo ist das einzelneurale Paradigma (mit seinen 
10 11 Einzelneuronen und ca. 10 15 Verbindungen zwischen diesen) ungeeignet. 
Infolgedessen wird der darauf basierte MED Ansatz durch die scheinbare 
Unerklarlichkeit beziehungsweise Zufalligkeit seiner klinischen Ergebnisse 

1 5 beeintrachtigt (fur transkraniale Magnetstimulation -TMS - wird auf den Oberblick 
[3] verwiesen). Weitere himaktivitatsverandernde MED Verfahren sind intra- and 
extrakraniale Elektroschocks. Alle MED Verfahren zielen darauf ab, 
Neuronmembranen zu depolarisieren, und solcherart vermehrt Aktionspotentiale 
zu erzeugen (sogenannte "Stimulation"). 

20 In CON-Terminologie lasst sich die elektromagnetische Depolarisation von 
Einzelneuronen durch MED-Verfahren als "Versklavung" bezeichnen (wobei 
sensorische Versklavung durch rhythmische Lichtblitze, akustische Signale, o.a. 
nicht betrachtet werden soil). Versklavung ist das Ergebnis der Anwendung 
kunstlicher Signale hoher Intensitat auf ein System, das dann wahrend dieser 

25 Anwendung das von auRen vorgegebene Signalmuster anstelle seines 
natiirlichen Verhaltens ubernimmt. In NA ist bekannt, dass unter schwacher 
Kopplung die Kommunikation in neuralen Netzen von der Kommensurabilitat der 
Frequenzen abhangt, mit denen die beteiligten Neuronen oszillieren. Demzufolge 
bedingt lokales Versklaven mit einer nicht oder nur schwach kommensurablen 

30 Frequenz die Ausschaltung der versklavten Neuronen aus ihrer normalen 
Kommunikation, und erzeugt solcherart eine virtuelle Lasion. Gegenwartige MED 
Verfahren betrachten virtuelle Lasionen als wesentlichen Effekt von TMS - der im 
einzelneuralen Paradigma unerklarlich ist. 

MED Verfahren betrachten die "Stimulation" als verbessert, wenn eine weniger 
35 hohe Intensitat benotigt wird. urn Einzelneuronen zum Feuern zu bewegen, 
beziehungsweise der von der Versklavung erfasste Raumbereich starker 
eingegrenzt werden kann. Wiederholte Versklavung am gleichen medizinisch 
sinnvollen Lokus kann den Anwendungszeitraum iiberdauernde physiologische 
Veranderungen bedingen, die in einigen Fallen positiv sind (z.B. statistisch. 
40 signifikante Verbesserung gewisser Typen klinischer Depressionen unter 
wiederholter TMS, siehe z.B. [3]). 
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In der Anwendung bei Menschen sind gegenwartige Verfahren nur in zweierlei 
Hinsicht individualisiert: erstens beim Auffinden eines geeigneten Lokus der zu 
depolarisierenden Einzelneuronen, zweitens beim Auffinden eines Bezugswertes 
fur die Starke des zur Versklavung fuhrenden elektromagnetischen Feldes ( 
5 angegeben in Prozent einer "motor threshold"). Bei einem neueren MED 
Verfahren [4] wird TMS diagnostisch eingesetzt, urn die individuelle Existenz oder 
Nichtexistenz physiologischer neuraler Verbindungen zwischen verschiedenen 
Hirngebieten, sowie, mit Einschrankungen, Ubertragungsgeschwindigkeiten 
zwischen verschiedenen Hirngebieten festzustellen. 

10 Gegenwartige MED Verfahren dienen vorzugsweise medizinischen 
Verwendungszwecken, in denen die Nachteile einer Versklavung (bis hin zum 
unwahrscheinlichen Fall von Krampfanfallen) durch deren medizinische Effekte 
zumindest kompensiert werden. Gegenwartige MED Verfahren umfassen 
lediglich passive Sicherheitsmaftnahmen wie Vermeidung hochfrequenter 

15 Versklavung, Ausschluss von Patienten mit Neigung zu epileptischen Anfallen, 
Abschalten des elektromagnetischen Feldes, falls ein Krampfanfall auftreten 
sollte. 

Zusammenfassend: Mit gegenwartigen Techniken und Verfahren kann eine 
20 modellbasierte, gesteuerte oder geregelte elektromagnetische Modifikation von 
Hirnaktivitat in vivo weder erreicht werden, noch bewegen sich diese Techniken 
und Verfahren auf diese Erreichbarkeit hin: MED ist im einzelneuralen Paradigma 
verhaftet, NA bewegt sich auf den Bau kunstlicher Systeme hin, und CON leidet 
an in-vivo Nichtanwendbarkeit bei komplexen heterogenen Systemen. 

25 

Aufgabe der vorliegenden Erfindung ist es, ein Verfahren und eine Vorrichtung fur 
die gesteuerte oder geregelte elektromagnetische Modifikation von Hirnaktivitat in 
vivo und die daraus resultierende Verhaltensmodifikation zu schaffen. 

30 Die Aufgabe wird durch Verfahren mit den Merkmalen im Anspruch 1 und durch 
eine Vorrichtung mit den Merkmalen im Anspruch 44 geldst Vorteilhafte 
Ausgestaltungen der Erfindung sind in den Unteranspruchen angegeben. 

Beim erfindungsgemaften Verfahren zur elektromagnetischen Modifikation von 
35 Hirnaktivitat wird eine Hirnaktivitatsmodell verwendet, das den Einfluss exogener 
elektrischer und/oder magnetischer Felder auf eine Hirnaktivitat beschreibt, und 
ein Verhaltensmodell verwendet, das den Zusammenhang zwischen Hirnaktivitat 
und Verhalten beschreibt. Hierdurch wird es erstmals moglich, gezielt mittels 
exogenen Inputs das Verhalten einer Person zu beeinflussen. 
40 Insbesondere betrifft die Erfindung die Verwendung eines Hirnaktivitatsmodells, 
sowie Ermittlung der individuellen, gegebenenfalls intraindividuellen und/oder 
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zeitabhangigen Nichtobservablen, sowie Ermittlung der individuellen, 
gegebenenfalls intraindividuellen und/oder zeitabhangigen 

Ubersetzungsoperatoren zwischen extrakranialem Signal und Kontrollkraft 
gestatten eine abgesicherte Intervention im Rahmen eines Steuer- oder 
5 Regelkreises, die wiederum ein zuverlassiges Erreichen eines Hirnaktivitatszieles 
bewirkt. Die Verwendung eines Verhaltensmodells gewahrleistet, dass das 
Erreichen eines Hirnaktivitatszieles das Erreichen des individuellen 
Verhaltensziels des Anwenders bewirkt. Der Erfindung liegt die Erkenntnis 
zugrunde, dass mit Verhaltensmodellen der Zusammenhang zwischen Verhalten 
10 und der Dynamik von HirnaktMtatscharakteristika quantifiziert wird, und mit 
geeignet individualisierten Hirnaktivitatsmodellen wiederum die Dynamik von 
Hirnaktivitatscharakteristika unter Angabe geeigneter Kontrollparameter 
quantifiziert wird, so dass eine zuverlassige Intervention zur Erreichung eines 
individuellen Verhaltenszieles moglich ist. 

15 

Die Erfindung wird nachfolgend beispielhaft anhand der Zeichnungen erlautert. In 
den Zeichnungen zeigen: 



Fig. 1 einen Transmitter in einer Schnittansicht, 
20 Fig. 2 den Transmitter aus Figur 1 in einer Ansicht von unten, 
Fig. 3 eine planare Projektion von Offnungen fur Sensoren und 

Transmitter entsprechend ihrer Anordnung an einem Helm, 
Fig. 4 schematisch einen Helm und Tragerachse nebst einer Kinnstiitze 

einer Kopfeinheit, und 

25 Fig. 5 eine weitere planare Projektion von Offnungen fur Sensoren und 
Transmitter entsprechend ihrer Anordnung an einem Helm, 
Fig. 6 prinzipielles Zusammenwirken von Kopfeinheit, Zwischeneinheit, 
Basiseinheit, 

Fig. 7 HAC-Modellmoden unter Variation von T(AuszCige), 
30 Fig. 8 HAC-Modellmoden unter Variation von a (Auszuge), 
Fig. 9 HAC-Modellmoden unter Variation von b (Auszuge), 
Fig. 1 0 HAC-Modellmoden unter Variation von c (Auszuge), 
Fig. 1 1 HAC-Modellmoden unter Variation von d (Auszuge), 
Fig. 12 ahnliche HAC-Modellmoden zu zwei verschiedenen Modellmoden I 
35 und II (Beispiel), 

Fig. 13 HAC-Modellmoden zu I unter sinusformigen exogenen Inputs 

verschiedener Starke (Beispiele), 
Fig. 14 HAC-Modellmoden zu II unter sinusformigen exogenen Inputs 
verschiedener Starke (Beispiele), 
40 Fig. 15 HAC-Modellmoden zu I unter exogenen Elementarloops niedriger 
Starke (Beispiele), 
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Fig. 16 HAC-Modellmoden zu II unter exogenen Elementarloops niedriger 
Starke (Beispiele), 

Fig. 17 HAC-Modellmoden zu I unter exogenen Elementarloops mittlerer 
Starke (Beispiele), 

5 Fig. 18 HAC-Modellmoden zu II unter exogenen Elementarloops mittlerer 
Starke (Beispiele), 

Fig. 19 eine Theta-HAC-Modellmode unter verschiedenen endogenen 

Inputs und einem exogenem Input (Beispiele), 
Fig. 20 Flussdiagramm des gesamten Verfahrens, 
1 0 Fig. 21 Flussdiagramm des Verfahrens S 2000 (Kalibrieren), 

Fig. 22 Flussdiagramm des Verfahrens S 2100 (Lokales Kalibrieren 1 
Teil), 

Fig. 23 Flussdiagramm des Verfahrens S 2100 (Lokales Kalibrieren 2 
Teil), 

1 5 Fig. 24 Flussdiagramm des Verfahrens S 2300 (Translokales Kalibrieren), 
Fig. 25 Flussdiagramm des Verfahrens S 3000 (Steuern und/oder Regeln 

von Hirnaktivitat), 
Fig. 26 ein Diagramm, das die Daten eines EEG-Kanals zeigt, 
Fig. 27 ein Diagramm, das einen Ausschnitt der Daten aus Fig. 26 im 
20 Phasenraum zeigt, 

Fig. 28 ein Diagramm, das im Phasenraum mit dem Wilson-Cowan Mod ell 

modellierte Neuronenaktivitat unter verschiedenen Inputs darstellt, 
Fig. 29 ein Diagramm, das im Phasenraum mit dem Wilson-Cowan-Modell 
modellierte Neuronenaktivitat bei gleichem Input unter 
25 verschiedenen Mischungswinkeln darstellt. 



Vorrichtung: 



Die erfindungsgemalie Vorrichtung 1 umfasst 

eine Oder mehrere Kopfeinheiten 2, verbunden mit 

einer Zwischeneinheit 3, von denen eine oder mehrere verbunden sind 

mit 

einer Basiseinheit 4 (Figur 6). 



35 Eine Kopfeinheit 2 umfasst ein Mellsystem mit Vorrichtungen zur 
elektromagnetischen Messdatengewinnung (..Sensoren"), Messdaten- 
vorverarbeitung, und Messdatenweitergabe an die Zwischeneinheit, sowie ein 
Stellsystem mit Vorrichtungen 5 zur Erzeugung von elektromagnetischen Feldern 
(..Transmitter"), sowie zur Umsetzung der von der Zwischeneinheit ausgehenden 

40 digitalen Steuerungs- bzw. Regelungsvorgaben in Transmittersignale (z. B. D/A- 
Umsetzereinheit). Eine Zwischeneinheit umfasst einen Rechner mit Software, mit 
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der die unten naher erlauterten Verfahren zur Steuerung von Hirnaktivitat bei 
nicht notwendiger Kalibrierung implementiert werden, sowie Verbindungen zu und 
von der Basiseinheit. Die Basiseinheit 4 umfasst einen Rechner, der erstens eine 
Datenbank m'rt den Modell- und Anwenderdaten enthalt, und auf dem zweitens 
5 die ubrigen unten naher dargestellten Verfahren implementiert sind. Die 
Zwischen- und die Basiseinheit sind auf dem gleichen oder auf verschiedenen 
Rechnern lokalisiert. 

Eine Ausgestaltung der Vorrichtung weist eine Kopfeinheit, eine Zwischeneinheit, 
10 und eine Basiseinheit auf. Die Kopfeinheit umfasst bspw. ein MeRsystem mit 
EEG-Kappe mit ihren extrakranialen Sensoren, Verbindungen zum Verstarker, 
Verstarker, Verbindungen zum A/D-Umsetzer, A/D-Umsetzer, Verbindungen zur 
Zwischeneinheit, sowie ein Stellsystem mit extrakranial befestigten 
stromfuhrenden Spulen als Transmittern, steuerbarer Stromversorgung fur diese 
15 Transmitter, nebst Verbindungen, D/A-Umsetzer, Verbindung von der 
Zwischeneinheit. Die Zwischeneinheit umfasst bspw. einen PC mit Bildschirm und 
Eingabetastatur, sowie Software, sowie Verbindungen zu und von der 
Basiseinheit. Die Basiseinheit umfasst bspw. eine leistungsfahige Rechnereinheit. 

20 Geeignete Sensoren sind beispielsweise EEG- oder MEG-Sensoren. Die MEG- 
Sensoren sind beispw. aus einem SQUID-Sensorelement mit geeigneter 
Auswertungseinrichtung zum Detektieren eines Magnetfeldes und 
Kiihleinrichtung ausgebildet. Die EEG-Sensoren weisen z. B. zwei Elektroden 
zum Messen einer elektrischen Potenzialdifferenz auf. 

25 

Eine vorteilhafte Ausgestaltung eines Sensors umfasst dessen teilweise oder 
vollstandige elektrische und/oder magnetische Abschirmung gegenuber seiner 
Umgebung, soweit dadurch dessen Funktion nicht behindert wird. 

30 Eine vorteilhafte Ausgestaltung der kopfnahen Teile des Meftsystems umfasst 
eine Vielzahl von Sensoren, die intra- und/oder extrakranial verteilt sind, diese 
Vielzahl von Sensoren wird als Sensorgitter bezeichnet. 

Eine vorteilhafte Ausgestaltung eines extrakranialen Sensorgitters umfasst 
35 dessen Fixierung beziiglich des Kraniums des jeweiligen Anwenders, sodass bei 
mehrfachem Auf- und Absetzen des Sensorgitters die Sensoren ihre jeweilige 
relative Position wieder einnehmen, beispielsweise durch Einpassen des 
Sensorgitters in einen Helm, dessen Innenseite die Kranialform des jeweiligen 
Anwenders nachbildet. 

40 
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Eine andere vorteilhafte Ausgestaltung des Sensorgitters umfasst implantierte 
Elektroden. 

Eine vorteilhafte Ausgestaltung eines extrakranialen Transmitters 5 umfasst eine 
stromfuhrende Spule 6 mit para-, dia-, oder ferromagnetischem Kern 7, wie in 
Figur 1 schematisch in Schnittansicht dargestellt, wobei die Pfeilrichtungen die 
Richtungen des Stromflusses symbolisieren. Der Transmitter 5 weist im 
Wesentlichen eine zylinderformige Form auf, wobei die Mantelflache und eine die 
Ruckseite bildende Stirnflache des Zylinders mit einer Abschirmung 8 
eingekleidet sind. An der von der Abschirmung freien Seite des Transmitters 
grenzen unmittelbar die Spule 6 und der Kern 7 an, und mit dieser Seite wird der 
Transmitter 5 im Betrieb auf das Kranium zur Abgabe exogener Magnetfelder 
ausgerichtet. An der Ruckseite des Transmitters 5 ist ein Halteelement 9 
angeoidnet, mit welchem der Transmitter 5 in einem Helm fixierbar ist. 
Der extrakraniale Transmitter 5 gegen Verformung geschiitzt sein, beispielsweise 
durch EingieRen der stromfuhrenden Telle in geeignetes Harz, oder Einbetten der 
stromfuhrenden Teile in stabiles Isoliermaterial. 

Der Transmitter 5 kann mit einer Kuhlvorrichtung versehen sein. 

Eine andere vorteilhafte Ausgestaltung eines intrakranialen Transmitters umfasst 
implantierte Elektroden. 

Eine vorteilhafte Ausgestaltung der kopfnahen Teile des Stellsystems umfasst 
eine Vielzahl von Transmittern, die intra- und/oder extrakranial verteilt sind, diese 
Anordnung von Transmittern wird als Transmittergitter bezeichnet. 

Eine vorteilhafte Ausgestaltung eines extrakranialen Transmittergitters umfasst 
dessen Fixierung bezuglich des Kraniums des jeweillgen Anwenders, sodass bei 
mehrfachem Auf- und Absetzen des Transmittergitters die Transmitter ihre 
jeweilige relative Position wieder einnehmen, beispielsweise durch Einpassen 
des Transmittergitters in einen Helm, dessen Innenseite die Kranialform des 
jeweiligen Anwenders nachbildet. 

In einer altemativen Ausgestaltung kann Fixierung von Sensor und/oder 
Transmittergitter auch kameraunterstutzt erfolgen, wobei die Position des Kopfes 
des Anwenders im Raum, sowie der Sensoren und/oder Transmitter bezuglich 
des Kopfes uber mehrere Kameras erfasst und real-time in 3D-Daten 
umgerechnetwird. 
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In einer weiteren Ausgestaltung werden als Sensoren und / Oder Transmitter 
intrakranial implantierte Elektroden verwendet, uber die sowohl EEG Messungen 
durchgefuhrt werden konnen, als auch Strome in das Hirn geleitet werden 
konnen. Zu diesen Elektroden fuhrende Leitungen und/oder deren Interfaces zur 
5 Rechnereinheit und/oder weitere Leitungen und/oder weitere Messgerate 
und/oder die zugehdrige Rechnereinheit und/oder der Energieversorger von 
Elektroden und/oder Rechnereinheit sind ebenfalls implantierbar, wodurch ein 
ambulanter Betrieb gestattet wird. 

Eine in Fig. 4 schematisch dargestellte vorteilhafte Ausgestaltung der kopfnahen 
10 Teile eines extrakranialen Mess- und Stellsystem umfasst einen auf seiner 
Innenseite die Kranialform des jeweiligen Anwenders nachbildenden Helm 10 mit 
durch eine Tragerachse 11 verlaufenden Verbindungskabeln und einer 
Kinnstutze 12. Sensor- und Transmittergitter im Inneren des Helms sind 
solcherart fixiert, dass sich beide Gitter Qberiappen - i.e. sich in der 
15 Nachbarschafl jedes Sensors hinreichend viele Transmitter befinden und 
umgekehrt. Eine planare Projektion der Uberiagerung des Transmitters mit dem 
Sensorgitter zeigt Figur 3 (hierbei sind Offnungen 13, Sensoren als Kreise und 
Offnungen 14 fur Transmitter 5 als Vierecke dargestellt). In der beschriebenen 
Ausgestaltung sitzt der Anwender auf einem Sessel mit Nackenstutze unterhalb 
20 des Helms 10. 

In einer altemativen Ausgestaltung ist das Sensorgitter intrakranial, und der Helm 
enthalt das extrakraniale Transmittergitter. 

25 In einer altemativen Ausgestaltung ist das Transmittergitter intrakranial, und der 
Helm enthalt das extrakraniale Sensorgitter. 

In einer altemativen Ausgestaltung sind sowohl Sensor-, als auch 
Transmittergitter intrakranial. 

30 

In einer vorteilhaften Ausgestaltung lasst sich die Sensordichte bzw. 
Sensorkonfiguration eines extrakranialen Sensorgitters einstellen. In einer 
weiteren vorteilhaften Ausgestaltung geschieht diese Veranderung automatisiert, 
gesteuert bzw. geregelt uber die Zwischeneinheit. 

35 

In einer vorteilhaften Ausgestaltung lassen sich die Transmitterdichte bzw. 
Transmitterkonfiguration eines extrakranialen Transmittergitters einstellen 
und/oder die Neigungswinkel jedes einzelnen Transmitters zum Kranium des 
Anwenders verandem. Eine planare Projektion einer mechanischen Halterung 
40 dieser Ausgestaltung zeigt Figur 5. Hierbei sind Offnungen fur Sensoren 13 als 
Kreise und Offnungen 14 fur Transmitter eckig dargestellt. Hier ist es moglich, 
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Transmitter 5 in den Offnungen 14 der Halterung zu verankern, und/oder 
Transmitter 5 gegenuber der Halterung zu kippen. Unter anderem konnen in 
dieser Ausgestaltung samtliche konventionellen Spulenkonfigurationen mit ihrer 
Anordnung, Ausrichtung und Feldrichtung dargestellt werden. 

Einrichtung, Funktionstests und Wartung der extrakranialen Vorrichtung erfolgen 
durch Techniker, die keine medizinische Ausbildung besitzen miissen. 

In einer vorteilhaften Ausgestaltung 1st die Vorrichtung mit konventionellem 
Schutz vor Stromausfallen und/oder Spannungsschwankungen versehen. 

Auf der Zwischeneinheit konnen real-time und automatisch Verfahren zur Dichte- 
und Positionierungsoptimierung von Sensoren und Transmittem, sowie 
zur Beseitigung von Artefakten ablaufen, das heiftt, nichtzufalligen 
Messwertverzerrungen, die durch kunstlich erzeugte Magnetfelder, 
Augenbewegungen, Muskelzuckungen u.a. entstehen. 



Verfahren: 

Das erfindungsgemalie Verfahren umfasst im einzelnen die nachstehend 
angegebenen, im folgenden naher erlauterten und mit „S" gekennzeichneten 
Verfahrensschritte. Ein- und Ausgangsdaten sind mit „D" gekennzeichnet. Die 
einzelnen Schritte werden im folgenden erlautert: 

Die Grundziige des erfindungsgemafien Verfahrens sind in einem 
Flussdiagramm 20 dargestellt, wobei Rechtecke Verfahrensschritte darstellen, 
Rauten Entscheidungsschritte, und Parallelogramme Input- oder Outputdaten von 
Verfahrensschritten. 

S 500 In diesem Schritt werden die fur ein zuverlassiges Erreichen des 
genannten Verhaltenszieles (D 2000) notwendigen Grolien sowie 
deren Dynamik mit Hilfe eines Verhaltensmodells (D 3000) 
spezifiziert. 

Das Verhaltensmodell stellt den Zusammenhang zwischen 
Verhalten und dem Verlauf bestimmter HAC dar. „HAC" steht im 
folgenden sowohl fur Singular als auch fur Plural von 
..Himaktivitatscharakteristikum". Die HAC werden aus Messsignalen 
nach bestimmten Vorschriften berechnet (D 2200). 
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Beispiele fur HAC sind von einem Sensor (Elektrode plus 
Referenzelektrode) gemessene Potentialdifferenzen, Oder ein 
Leistungsspektrum des vom Sensor innerhalb eines Zeitfensters 
empfangenen Signals, oder der Quotient von Beta (13-30 Hz) - 
EEG-Aktivitat zu Alpha (8-12 Hz) -EEG-Aktivitat fur diesen Sensor 
(s.a. [11]). 

Beispiel fur ein Verhaltensmodellsind das reduzierte Davidson- 
Modell (siehe z.B. [11]), wonach positive Emotionen mit einem 
hoheren Quotienten von Beta (13-30 Hz) - EEG-Aktivitat zu Alpha 
(8-12 Hz) - EEG-Aktivitat im linksfrontalen Kortex einhergehen. 

Das Verhaltensmodell (D 3000) gestattet es, das Verhaltensziel (D 
2000), welches grundsatzlich die Dauemaftigkeit des gewiinschten 
Verhaltens umfasst, in einen Soll-Verlauf des oder der HAC (D 
2100) umzuwandeln, und somit zu quantifizieren. 

Aus Art und Anzahl der HAC ergeben sich Mindestanforderungen 
an die Apparatur (D 4000). Beispielsweise verwendet man zum 
Errechnen des Quotienten von Beta (13-30 Hz) - EEG-Aktivitat zu 
Alpha (8-12 Hz) - EEG-Aktivitat im linksfrontalen Kortex zumindest 
eine EEG-Elektrode plus eine Referenzelektrode. Zum kontrollierten 
Erreichen des Soll-Verlaufs des o.g. HAC verwendet man zumindest 
einen Transmitter. 

Die verwendeten HAC sind aus den Variablen des verwendeten 
generischen Himaktivitatsmodells (D 1000) zu berechnen. 

Beispiel: Auf Fourier-Zerlegung basierende HAC fur ein umgrenztes 
Himgebiet sind aus der Dynamik der neuralen Oszillatoren des 
Wilson-Cowan-Modells (siehe z.B. [1]) berechenbar. 

Aus den HAC ergeben sich Mindest-Kalibrierungsanforderungen (D 
1100) bezuglich des verwendeten generischen Modells. Je weniger 
robust sich die Berechnung der HAC gestaltet, desto starkere 
Anforderungen sind an die Kalibrierung zu stellen. Falls die HAC auf 
gekoppelten, nicht-robusten GroBen verschiedener ortlicher 
Herkunft basieren, ist die lokale Kalibrierung durch translokale zu 
erganzen. In jedem Fall sind Art (bspw. „keine, lokal, translokal", D 
1110 in Fig. 21) und Umfang (bspw. Liste der Sensoren sowie 
Transmitter, D 1120) zu benennen (E 100). 



20 



25 
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Beispiele: Fur den robusten, nichtlinearen Ahnlichkeitsindex (siehe 
bspw. [5]) ist eine Kalibrierung nicht unbedingt erforderlich. Fur den 
o.g. Beta/Alpha-Quotienten ist eine lokale Kalibrierung des Wilson- 
5 Cowan-Modells erforderlich. Bei komplexen kognitiven Leistungen 

(bspw. musischer Kreativitat, s. [19]) ist in der Regel translokale 
Kalibrierung erforderlich. 

S2000 Kalibrieren umfasst das Ermitteln individueller, gegebenenfalls 
10 zeitspezifischer Werte fur die im generischen Himaktivitatsmodell 

aus Populationsmittelwerten geschatzten Parameter, endogenen 
Input, sowie die Quantifizierung des ebenfalls individuellen, 
gegebenenfalls zeitspezifischen, exogenen Einflusses 
elektromagnetischer Felder auf Hirnaktivitat. („Elektromagnetisch" 
1 5 steht fur ..elektrisch und/oder magnetisch" und wird im folgenden mit 

„EM" abgekurzt.) Umfang und Qualitat der Kalibrierung beeinflussen 
die Zuverlassigkeit der Modifikation ein oder mehrerer 
Himaktivitaten. 

20 Aus dem Kalibrieren resultiert ein spezifisches Hirnmodell (D 1200), 

das exogenen Input als Kontrollvariable innerhalb des Modells 
quantitativ erfasst (D 1300). Kalibrieren ist naher ausgefuhrt in Figur 
21. 

25 Man unterscheidet im einfachsten Fall lokales (siehe Figuren 22-23) 

und translokales (siehe Figur 24) Kalibrieren. 

Erstes Ziel der lokalen Kalibrierung ist es, fur jeden Sensor bei 
einem bestimmten Anwender in einem bestimmten Zeitintervall 
30 einzugrenzen, welche Himaktivitaten im Erfassungsbereich des 

Sensors stattfinden (Ablauf dargestellt in Figur 22). 

Zweites Ziel der lokalen Kalibrierung ist es, fur jede Hirnaktivitat im 
Messbereich jedes Sensors herauszufinden, wie diese durch von 
35 zumindest einem Transmitter erzeugte exogene EM Felder 

beeinflusstwerden kann (Ablauf dargestellt in Figur 23). 

Ziel translokaler Kalibrierung (S 2400) ist es, Art und Umfang des 
Einflusses der Aktivitat anderer Hirnbestandteile auf die Aktivitat des 
40 betrachteten Himbestandteils festzustellen. 
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Das Kalibrieren wird nachfolgend anhand der Fig. 21 bis 24 erlautert: 
Fig. 21 zeigt einen Ausschnitt des Verfahrens, bei dem mit den bereits oben 
erlauterten Daten D 1100, D1110 und D 1120 die Entscheidung gemalX dem 
oben beschriebenen Entscheidungsschrittes E 100 getroffen wird, ob eine 
5 Kalibrierung durchgefuhrt werden soli oder nicht. Soil keine Kalibrierung 
durchgefiihrt werden, gent das Verfahren auf den Schritt S 3000 uber, wie es 
oben bereits ausgefuhrt ist. Soil hingegen eine Kalibrierung durchgefuhrt werden, 
geht das Verfahren auf den Schritt S 2100 uber, bei welchem ein lokales 
Kalibrieren ausgefuhrt wird, das unten naher anhand der Fig. 22 bis 23 erortert 
10 wird. 

Da beim lokalen Kalibrieren Himaktivitaten lokal kalibriert werden, kann es sein, 
dass ein und dieselbe Hirnaktivitat mehrfach von unterschied lichen Sensoren 
bzw. Sensoreinheiten erfasst und kalibriert werden. Deshalb ist es zweckmaliig, 
15 im Schritt S 2300 identische Himaktivitaten als solche zu identifizieren, wobei 
mehrfach kalibrierte Himaktivitaten nur einmal beriicksichtigt werden oder 
mehrere Kalibrierungen einer Hirnaktivitat kombiniert ausgewertet werden. Die 
entsprechenden Daten werden dann in die Datensatze D 1200 und D 1300 
integriert. 

20 

Falls eine translokale Kalibrierung auf D 1110 erforderlich sein sollte, wird dies im 
Schritt E 120 gepruft und gegebenenfalls eine translokale Kalibrierung S 2400 
ausgefuhrt, die unten naher anhand von Fig.24 erlautert wird. Die 
entsprechenden Daten werden dann in die Datensatze D 1200 und 1300 
25 eingefiigt. 

Die lokale Kalibrierung wird nachfolgend anhand von Fig. 22 und 23 erlautert. 

Hierzu wird mit jedem Sensor EM Hirnaktivitat gemessen (S 2110). Die 
30 entstehende Zeitreihe (D 2205) wird in einen Ist-Verlauf von HAC umgerechnet 
(S 2120), wodurch sich Ist-Verlaufe von HAC ergeben (D 2210). 

Ein Hirnaktivitatsmodeil (D 1000) wird in der Regel mittels Differentialgleichungen 
beschrieben, in die exogener Input unter bestimmten Einflussannahmen (D 1049) 
35 integriert wird. 

Eine Modellmode ist eine Losung der Differentialgleichungen des solcherart 
erweiterten Hirnaktivitatsmodells fur einen Satz von Nichtobservablen, d.h. von 
Parametem, endogenen und exogenen Inputs. 

40 
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Sie werden im Schritt S 2115 berechnet und in einer Datenbank (D 1051) 
abgespeichert. Diese Datenbank enthalt die Zuordnung von Satzen von 
nichtobservablen Modellmoden. Mit Hilfe der in der Datenbank D 1051 
gespeicherten Daten und/oder unmittelbar im Schritt S 2115 unmittelbar 
5 berechneten Daten werden die HAC-Verlaufe der Modellmoden berechnet (S 
2120). Hieraus ergeben sich theoretisch mogliche HAC-Verlaufe (D 1030). 

Da HAC i.d.R. vereinfachend wirken, ist es moglich, dass mehrere Modellmoden 
in ahnlichen (d.h. innerhalb von Fehlerschranken nicht unterscheidbaren) HAC- 

10 Modellmoden resultieren. Ebenso ist es moglich, dass mehrere Satze von 
Nichtobservablen in der gleichen Modellmode resultieren. Dieser Zusammenhang 
zwischen HAOModellmoden, Modellmoden, und Satzen von Nichtobservablen 
tritt analog auch in Eichtheorien auf, wie bspw. Elektrodynamik (entsprechend 
einem HAC hier bspw. „Energiedichte in einem Raumsektor", entsprechend der 

15 Modellmoden verschiedene elektrische und magnetische Felder, die diese 
Energiedichte generieren, und entsprechend den Satzen von Nichtobservablen 
verschiedene Potentiate, die wiederum den elektrischen/magnetischen Feldem 
zu Grunde liegen konnen). Fur stochastische Moden sind beispielsweise deren p- 
te Momente im Zeitverlauf gespeichert Modellmoden und ggf. HAC-Moden sind 

20 in einer Datenbankgespeichert und/oder schnell berechenbar. Es sind nur endlich 
viele Satze von Nichtobservablen gespeichert (Gitter), desgleichen fur 
Modellmoden und HAC-Modellmoden. Bei Bedarf besteht die Moglichkeit, das 
Gitter zu verfeinern und/oder zu erweitern (bspw. durch neue Parametersatze 
und/oder neue Arten und/oder neue Koeffizienten endogenen und/oder exogenen 

25 Inputs). 

Es ist moglich, dass eine Mode auf dem Weg von ihrem Trager zum jeweiligen 
Sensor Abschwachungen und/oder Verzerrungen und/oder 
Phasenverschiebungen unteriiegt. (Als „Trager* einer Mode wird dasjenige 
30 Neuronenensemble definiert, dessen Aktivitat die betreffende Mode erzeugt). 
Das am jeweiligen Sensor ankommende Bild der Mode wird als Bildmode 
bezeichnet, analog fur zugehorige HAC-Moden. 

In dem o.g. Ist-Verlauf von HAC werden die gemali a-priori Annahmen (D 1040, 
Figur 22) gegebenenfalls veranderten Bilder von HAC-Modellmoden identifiziert 

35 (S 2130). Im einfachsten Fall ist das Bild einer HAC-Modellmode proportional zur 
HAC-Modellmode selbst und kann (aufler in Fallen vernachlassigbarer 
Amplitude) mit dem Quotienten aus HAC-Modellmode und dem Reziproken eines 
Abschwachungsfaktors identifiziert werden. Dieser einfachste Fall wird im 
folgenden zugrundegelegt. Eine (HAC-) Modellmode ist im einfachsten Fall lokal 

40 und beschreibt damit eine lokalisierte Himaktivitat. 
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Dekomposition gemafi S 2130 (s. Figur 22) beinhaltet das Auffinden von HAC- 
Modellmoden (D 1030) im empirischen HAC-Verlauf, das Identifizieren von 
diesen HAC-Modellmoden moglicherweise zugrundeliegenden Modellmoden und 
das Identifizieren von diesen Modellmoden moglicherweise zugrundeliegenden 
5 Satzen von Nichtobservablen, wobei hier nur Satze von Nichtobservablen mit 
exogenem Input gleich Null zugelassen sind. Endogene Inputs werden im 
einfachsten Fall fur die Zwecke des lokalen Kalibrierens als konstant 
angenommen, im zweiteinfachsten Fall als sinusformig, im dritteinfachsten Fall 
als oszillierend. 

10 Hilfreich fur die Praktikabilitat des Verfahrens ist die Beschrankung auf die 
Untersuchung relevanter HAC-Modellmoden. Eine in den Ist-Verlaufen 
aufgefundene HAC-Modellmode soil relevant" heifien, wenn sie im Ist-HAC- 
Verlauf oberhalb vordefinierter Schwellenwerte fur Rauschen, Messfehler, o.a. 
festgestellt werden kann. 

15 

Ergebnis der Dekomposition sind eine Liste relevanter HAC-Modellmoden (D 
2300) im Messbereich des jeweiligen Sensors, eine Liste moglicher Modellmoden 
zu jedem Element von D 2300 (D 2350) und eine Liste moglicher 
Parametersatze, sowie mogliche endogene Inputs zu jeder Modellmode in D 
20 2350 (D 2400). 

Wie in Fig. 23 dargestellt, werden fur jeden Sensor Si und fur jeden Transmitter Tk 
die folgenden Verfahrensschritte ausgefuhrt: 

25 Aus den Parametersatzen mit endogenem Input (D 2400) werden nun durch 
Testen (S 2170, Senden von Testsignalen mit Hilfe von Transmittern) diejenigen 
Satze (D 2410) von Nichtobservablen herausgefunden, die die jeweilige HAC- 
Modellmode mit Testsignal erklaren. Dieses umfasst es, mit Hilfe des Testsignals 
zu ermitteln, wie sich das Transmittersignal in exogenen Input ubersetzt, das 

30 heilit, in einen Bestandteil der Modellgleichungen. Sobald dieses geklart ist, wird 
der Einfluss des exogenen Signals auf die entsprechende Modellmode respektive 
HAC-Modellmode berechenbar und zur tatsachlichen Modifikation von 
Himaktivitat (S 3000) verwendet. 

35 Ein Testsignal (D 2500) ergibt sich im einfachsten Fall aus den moglichen 
endogenen und exogenen , Inputs (umfassend: konstant, oszillierend, mit 
Elementarloops, dh. zeitverzogerter Verwendung eines Messsignals als 
Transmittersignal, und viele andere). 

40 Das Testsignal wird in das generische Hirnmodell wie folgt integriert: man ersetze 
„lnput u in den Modellgleichungen durch „endogenen plus exogenen Input". 
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Exogener Input := Ubersetzungsoperator(Transmittersignal) 

Der Ubersetzungsoperator 0 stellt den Zusammenhang zwischen exogenem, 
5 physikalisch messbarem Signal und Kontrollvariablen in den Gleichungen des 
spezifischen, d. h. kalibrierten, Hirnmodells dar. Je nach Parametersatz und 
endogenem Input, die der einzelnen Mode unterliegen, sowie Gestalt und 
Amplitude des Transmittersignals, verandern verschiedene Werte der 
Kontrollvariable die betreffende Mode, was wiederum physikalisch messbar ist. 

10 

Ein Testsignal hat die Aufgabe, zwischen verschiedenen Satzen von Parametern, 
endogenen Inputs, und ggf. Gbersetzungsoperatoren fur die betrachtete HAC- 
Mode zu trennen. Ein Testsignal ist (wie jedes Transmitlersignal) gekennzeichnet 
durch eine exteme Referenzgrolie, sowie eine i.d.R. zeitabhangige Gestalt am 
15 Ort des Transmitters, und eine - moglicherweise einelementige - Folge von 
Amplitudenmultiplikatoren. 

Jedes Transmittersignal und damit auch jedes Testsignal wird hier i.d.R. in 
Hybridnotation angegeben, d.h. ausgedruckt in den Variablen des 
20 Hirnaktivitatsmodells, mit einer Folge von Amplitudenmultiplikatoren, sowie unter 
Angabe einer externen, die Signalerzeugung spezifizierenden Kenngrode. 

Beispielsweise kann das extrakraniale Magnetfeld mit 0.01 Tesla mit Hilfe einer 
von einem Strom der Starke l 0 durchflossenen Spule erzeugt werden. Wenn die 

25 Gestalt der bezuglich einer gewahlten Spannungseinheit auf Amplitude 1 
normierten Gestalt der durch die zeitlichen Anderungen des Magnetfelds 
induzierten Spannung sin(2*Pi*t*5/1000) ist (die Spannung ist als Variable des 
Hirnaktivitatsmodells angenommen, t ist die Zeit in Millisekunden), so ergibt sich 
bei sonst gleichen Bedingungen bei Stromen der Starke lo, 2*l 0 , 4*l 0 , 8*l 0 

30 folgende Signalnotation: 

(0.01T, sin(2*Pi*t*5/1000), (1,2,4,8)). 

Weiterhin kann mit z(t) als auf Spannungsamplitude 1 normiertem 
Systemverhalten sich ein urn 10 Millisekunden ruckgekoppeltes Signal mit 
35 Stromstarken analog zu vorher ergeben als:(0.01T, z(t-10), (1 ,2,4,8)). 

Die Zeitdauer der Teilsignale (bei l 0 , bei zwei mal lo, etc.) richtet sich nach dem 
gewunschten Effekt (bspw. transient versus Grenzzyklus), ebenso die moglichen 
Pausen zwischen je zwei Teilsignalen. 

40 
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Der Ubergang vom Transmittersignal zum Signal am Trager der jeweiligen Mode 
ist Gegenstand von Einflussinvarianten und -annahmen (D 1050). 

Die moglichen Amplitudensequenzen unterliegen technischen und/oder 
5 gesundheitlichen Rahmenbedingungen (D 4100, D 4200). 

Im einfachsten Fall gibt es beispielsweise zwei verschiedene Satze von 
Parametern und endogenen Inputs, die zu verschiedenen Modellmoden fiihren, 
aber zu ahnlichen HAC-Modellmoden. Das heifit, durch blofte Betrachtung von 

10 HAC kann nicht entschieden werden, welcher Satz von Nichtobservablen der 
betreffenden Hirnaktivitat unterliegt Also wird zu einer Gestalt des Testsignals 
eine Folge von Amplitudenmultiplikatoren rnit Hilfe des generischen Hirnmodells 
durch Einsetzen der beiden Satze numerisch solcherart berechnet (S 2160), dass 
unter a priori Annahmen fur den Ubersetzungsoperator (bspw. Kombination von 

15 Phasenverschiebung und Abschwachung des Signals) mit Hilfe des Testsignals 
Modellmoden erzeugt werden, deren resultierende HAC-Modellmoden fur beide 
Satze verschieden sind. 

Aus dieser numerischen Berechnung ergibt sich i.d.R. eine Menge mdglicher 
20 Gestalten fur Testsignale, mit jeweils einer Folge moglicher 
Amplitudenmultiplikatoren. Aus dieser Menge wird dann mit Hilfe vordefinierter 
Kriterien (bspw. minimale Storung aller auftretenden Hirnaktivitaten mit 
Ausnahme der jeweils zu untersuchenden Mode durch das Testsignal) ein 
optimales Testsignal ausgewahlt. 

25 

Man beginnt mit dem Senden des Signals (S 2170) mit der ersten Amplitude der 
Sequenz, misst mindestens mit dem betrachteten Sensor S\ und hochstens mit 
alien Sensoren des Sensorgitters (s. z. B. Figur 5) die elektromagnetische 
Hirnaktivitat (S 21 10), erhalt daraus Zeitreihen (D 2520), aus denen man Veriaufe 
30 von HAC errechnet (S 21 20). 

Anschliefiend erfolgt die Dekomposition des Ist-HAC-Verlaufs nach HAC- 
Modellmoden mit exogenem Input (S 2135). Hierfur sind nur Parametersatze und 
endogene Inputs aus D 2400 (s. Figur 22) zugelassen, wobei nun an Stelle des 
35 dortigen „exogener Input gieich Null" ein moglicherweise anderer exogener Input 
tritt, dessen Gestalt sich aus den Einflussinvarianten und -annahmen ergibt, und 
dessen Starke sich i.d.R. aus den (i.d.R. nichtlinearen) Veranderungen bei 
verschiedenen Amplituden einer Amplitudensequenz ergibt. 

40 Falls diese Dekomposition mehrere bestapproximierende Satze von 
Nichtobservablen liefert (E 210), so ist ein weiteres Testsignal zu errechnen, und 



WO 03/085546 



-20- 



PC17EP03/03545 



die auf S 2160 folgende Schleife erneut zu durchlaufen, wobei dann bei S 2135 
nicht mehr die Satze von Parametem und endogenen Inputs aus D 2400 
beriicksichtigt werden, sondern D 2410 (s. Figur 23). Liefert mehrfaches 
Durchlaufen der Schleife keine Ergebnisse, so ist die Ausrichtung des 
5 Transmitters zu verandern Oder ein anderer Transmitter zu wahlen. 

Insgesamt ist nach Abschluss der lokalen Kalibrierung bekannt, 

- welche Modellmoden im Messbereich eines Sensors liegen, 

- welcher Parametersatz und welcher endogene Input die 
1 0 jeweilige Modellmode erzeugen, 

- wie die jeweilige Modellmode zu beeinflussen ist. 

Zusatzlich ist es im Anschluss an lokales Kalibrieren im Schritt S 2300 moglich, 
identische Trager von Modellmoden zu identifizieren. 

15 Fur Ari als eine dem Sensor Sk zugeordnete Mode, gekennzeichnet durch 
Parameter P w , endogenen Input lw(t) und UbersetzungsoperatorenC7 Wi;r (t), und 

Amn(t) als eine dem Sensor S m zugeordnete Mode, gekennzeichnet durch 
Parameter P mn , endogenen Input l mn (t) und Ubersetzungsoperatoren U mnx (t) (% 

Nummerierung der Transmitter) gilt: wenn alle P w gleich den entsprechenden P mn 
20 sind, und l w (t) gleich l mn (t) ,dann sind beide Moden identisch. Diese identischen 
Moden haben nicht denselben Trager, wenn es einen Transmitter T x gibt, so 

dass &„ a (\) ungleich U^X) und/oder U mix {\) ungleich U mn%x (X) ist. 

Verfahrensschritt S 2300 ist nicht essentiell, bewirkt aber eine mogliche 
Verkleinerung der Funktionalmatrix (D 2410, s. Figur 24) beim im folgenden 
25 dargestellten translokalen Kalibrieren. 

Mit der translokalen Kalibrierung (S 2400) werden Art und Umfang des Einflusses 
anderer Moden und/oder sensorischen Inputs auf eine Mode festgestellt 
(Indirekte Fremdeinwirkung, i.e. Darbietung visueller, akustischer, taktiler, und 

30 anderer Reize erfolgt uber das sensorische System, bedingt Verzogerungen, und 
wird als Teil des endogenen Inputs modelliert). Werkzeug hierfur ist die 
„Funktionalmatrix" (D2410): wenn insgesamt n verschiedene Modellmoden im 
Messbereich mindestens eines der Sensoren aufgefunden wurden, so ist diese 
Funktionalmatrix eine n*n Matrix , in deren (ij)-ter Zelle ein Eintrag zu finden ist, 

35 wenn die Mode i auf die Mode j einen nachweisbaren Einfluss hat. Der Eintrag 
quantifiziert diesen Einfluss, sodass in den Gleichungen des Hirnmodells fur die 
Mode j „endogener Input" ersetzt wird durch „gegebenenfalls zeitverzogerte 
Funktion der i-ten Mode plus sonstiger endogener Input". 

40 Das Verfahren zum Fullen benotigter Teile der Funktionalmatrix ist in Figur 24 
dargestellt. Kenntnis der Funktionalmatrix oder von Teilen der Funktionalmatrix 
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gestattet die Berechnung von mehre Konstltuenten des betreffenden 
Himaktivitatsmodell umfassenden Phanomenen (bspw. Synchronisierung, Phase- 
Locking, und vieles mehr), wobei konventionelle numerische Verfahren, neurale 
Netze, nichtlineare Dynamik, und ahnliches verwendet werden. 

5 

Die aus dem lokalen Kalibrieren erhaltenen Ubersetzungsoperatoren (D 1300) 
sowie weitere Nichtobservable gestatten fur alle Moden eine HAC-Prognose (S 
2410), die translokale Einflusse nur summarisch berucksichtigt. Bei aktivem 
Andern einer Mode (S 2420) wird ein moglicher Unterschied zwischen Ist-HAC 
1 0 und Prognose-HAC feststellbar. Diese wird im Schritt S 2430 erfasst. Diese 

Unterschiede beruhen im Wesentlichen auf translokalen Einflussen. Im Schritt S 
2440 werden diejenigen Funktionalmatrizen berechnet, deren Einbeziehung in die 
Gleichungen des Hirnaktivitatsmodells geeignet sind, den erfassten Unterschied 
zu erklaren. 

15 

In E 400 wird gepruft, ob diese Zuordnung von Unterschied auf Funktionalmatrix 
eindeutig ist. 1st dieses nicht der Fall, so erfolgt in Schritt S 2221 eine 
Modifizierung des Testsignals zurweiteren Einengung der Menge moglicher 
Funktionalmatrizen. Diese Iteration wird abgebrochen, sobald zumindest die 
20 gewunschten Elemente der Funktionalmatrix eindeutig feststehen. 

Im Anschluss an die Kalibrierung hat das vorzugsweise EM Steuem oder Regeln 
(S 3000) von Hirnaktivitat uber himaktivitatsmodellbasierte Riickkopplung hat das 
Erreichen sowie Beibehalten von Soll-Verlaufen von HAC zum Ziel. Dies ist in 

25 Fig. 25 im Detail dargestellt. 

Hierbei werden zunachst Transmittersignale berechnet, die theoretisch geeignet 
sind, den Ist-Verlauf der HAC (D 1010) in den Soll-Verlauf zu uberfuhren (S 
3100). Diese Transmittersignale setzen sich vorzugsweise aus einfachen 
Signalen mit im Rahmen des spezifischen Hirnmodells berechenbarer und ggf. im 

30 Rahmen der Kalibrierung bereits getesteter Wirkung zusammen, wobei sowohl 
Mehrfachwirkungen (Auswirkungen von Signalen eines Transmitters auf mehrere 
Hirnaktivitaten), als auch raumliche Zusammensetzung (Signale mehrerer 
Transmitter), als auch zeitliche Zusammensetzung (Sequenz von Signalen) in 
Frage kommen. Nach Machbarkeitsprufung der errechneten Signale (bezuglich 

35 Gesundheitsgrenzwerten (D 4200), Apparaturbeschrankungen (D 4100), ggf. 
Ausbreitungseinschrankungen (D 4300) wird mit Hirfe einer Nutzenfunktion (z.B.: 
kleinstmogliche Feldstarke pro Transmitter) das beste Transmittersignal 
ausgewahlt und gesendet (S 3400). Die theoretische Wirkung des Signals ist in 
einem spezifischen Himaktivitatsmodell mit Einflusskalibrierung bekannt und 

40 gestattet HAC-Verlaufsprognosen (D 4040). Falls im Rahmen der laufenden EM 
Messung und Ist-HAC-Berechnung bei einem Vergleich von Ist- und Soll-HAC- 
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Verlaufen (S 3500) deutliche Abweichungen auftreten (E 300), so wird S 3000 
vorlaufig unterbrochen, urn die Kalibrierung mit einem weiteren 
Kalibrierungsvorgang (S 2000) zu verbessem. Anschlieliend geht das Verfahren 
wieder auf den Schritt S 3100 uber, bei dem die Transmittersignale neu 
5 berechnet werden. 

Die wesentlichen Verfahrensschritte werden nun anhand einiger Beispiele naher 
ausgefuhrt. Die Beispiele werden anhand von Ablaufen durch die 
Flussdiagramme (Figur 20 - 25) erlautert. 

10 

Beispiel 1: 

Hour 20: 

1 5 D 2000 (Verhaltensziel): 

Verstarken positiver Emotionen fur die Dauer der Anwendung. 

D 3000 (Verhaltensmodell): 

Als reduziertes Davidson-Modell (siehe z.B. [11]), wird foigender Zusammenhang 
20 bezeichnet: positive Emotionen gehen mit einem hoheren Quotienten von Beta 
(13-30 Hz) - EEG-Aktivitat zu Alpha (8-12 Hz) - EEG-Aktivitat im linksfrontalen 
Kortex einher. 

S 500 (Spezifizieren HAC und deren Berechnung. Spezifizieren Soll- 
25 Verlaufe von HAC): 

Das Davidson-Modell basiert auf dem Leistungsspektrum von EEG-Signalen. Ein 
geeignetes HAC ist somit eine Folge quadrierter Absolutbetrage von 
Fourierkoeffizienten, die uber Fast Fourier Transformation berechnet wird, 
beispielsweise zu den Frequenzen von 1 bis 50 Hertz. Abgeleitetes HAC ist 

30 12 

30 pos\^Yafi' f^Lft 1 • (Annahme: nur ganzzahlige Frequenzen, Fourierfenster 

beispielsweise 500 Millisekunden, Absolutbetrag des Fourierkoeffizienten der i- 
Hertz-Mode. Urn die Frequenzen von 1 bis 50 Hz zu uberdecken, sind nach 
Nyquist mindestens 100 Messwerte erforderlich). HAC(100) ist das 
Leistungsspektrum bezuglich der Messwerte 1 bis 100, HAC(101) das 
35 Leistungsspektrum bezuglich der Messwerte 2 bis 101, etc. Pos(100) ist aus 
HAC(100) abgeleitet, Pos(101) aus HAC(101) etc. Dieses abgeleitete HAC soil 
erhoht werden, also beispielsweise: ab einem Zeitpunkt bzw. dem 
entsprechenden Messpunkt soil gelten: Pos(t, mit Einfluss) > 2*Pos(t, ohne 
Einfluss) . Diese Anforderung bestimmt alle Soll-Verlaufe des HAC. 

40 
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D 4000 (Minimale Apparaturvoraussetzungen): 

Oberflachen-EEG mit einer Elektrode linksfrontal, einer Referenzelektrode z.B. 
am Ohr (s. [16]), einer extrakranialen Spule (als Transmitter) linksfrontal. Mit 100 
Messwerten pro halber Sekunde (s. S 500) ergibt sich eine minimale Sampling- 
5 Rate von 200/sek. 

D 2100 (Soll-Verlauf von HAC): 

Pos(t, mit Einfluss) > 2*Pos(t, ohne Einfluss) 

1 0 D 2200 (HAC-Berechnungsvorschriften): 

Mit Hilfe der Fast Fourier Transformation werden die quadrierten Absolutbetrage 
fi 2 der Fourierkoeffizienten der i-ten Mode (i zwischen 1 und 50 Hertz) berechnet. 

30 12 
1=13 /=8 

1 5 D 1000 (Generisches Hirnaktivitatsmodell): 

Vereinfachtes Wilson-Cowan-Modell (abgeleitet aus [1], wobei Refraktarzeiten 
auf Null gesetzt wurden, ebenso Leitungsverzogerungen, u.v.m.), wobei 
Hirnaktivitat auf der Aktivitat neuraler Oszillatoren basiert. Ein neuraler Oszillator 
oszilliert oder oszilliert nicht, abhangig von seinem Input, sowie ebenfalls nicht 

20 observablen physiologischen Parametem. Ein neuraler Oszillator ist ein 
Ensemble aus miteinander verbundenen erregenden und hemmenden 
Einzelneuronen (in [1] von der Grolienordnung 10 5 Einzelneuronen, wobei hier 
auch Vielfache davon zugelassen sind, sofern das jeweilige Neuronenensemble 
durch die folgenden Gleichungen beschreibbar ist), modelliert durch ein System 

25 zweier nichtlinearer Differentialgleichungen: 

rx = -x+ S(ax-by +P X ) 
ry = -y + S(cx-dy + p y ) . 

Hierbei ist #0?) :=l/(l + exp(5-£)) eine sigmoide Funktion. „x", sind 
30 elektromagnetische GroRen, die in der Form otx + py observabel sind, wobei a, p 
Mischungswinkel heifien und aus dem Interval! [-1,1] stammen, „x" fur die 
Gesamtheit der erregenden Einzelneuronen, „y" fur die der hemmenden 
(Gesamtheit = Gesamtheit innerhalb des neuralen Oszillators). „a", „b tt , „c M f „d M , 
sind langsame Nichtobservable. „a" stellt den Einfluss der Gesamtheit der 
35 erregenden Einzelneuronen auf sich selbst dar, analog „d" fur die Gesamtheit der 
hemmenden Einzelneuronen. „b" stellt den Einfluss der Gesamtheit der 
hemmenden Einzelneuronen auf die Gesamtheit der erregenden Einzelneuronen 
dar, ,,0" stellt den Einfluss der Gesamtheit der erregenden Einzelneuronen auf die 
Gesamtheit der hemmenden Einzelneuronen dar. n p x u stellt den endogenen 

40 Input von aufterhalb des neuralen Oszillators in die Gesamtheit der erregenden 
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Einzelneuronen dar. „p y " stellt den endogenen Input von aufcerhalb des 

neuralen Oszillators in die Gesamtheit der hemmenden Einzelneuronen dar. 
Inputvariable werden als innerhalb des Modells nicht notwendigerweise als 
langsam angenommen, auf diesen Aspekt wird bei der Erlauterung translokaler 
5 Kalibrierung genau eingegangen. Alle Zeiten werden in Millisekunden 
ausgedruckt, alle anderen GroBen sind reine Zahlen. Das r auf der linken Seite 
der vereinfachten Wilson-Cowan-Gleichungen gibt die Membranzeitkonstante in 
Millisekunden an, eine weitere Nichtobservable. 

Urn Fremdeinwirkung zu umfassen, wird das vereinfachte Wilson-Cowan-Modell 
10 erweitert: Direkte elektromagnetische Fremdeinwirkung wird als Teil des Inputs 
verstanden: p x = p Xttndo + p x%exo und p y = p yt6ndo + , wobei 

p^^u^Transmittersignal) den exogenen Input in den betrachteten neuralen 
Oszillator auf Grund des vom i-ten Transmitter ausgestrahlten Signals 
bezeichnet, und f ,endo M fur ..endogen" steht U ist ein Ubersetzungsoperator, der 
15 unter anderem von der Entfernung des betrachteten neuralen Oszillators vom 
betreffenden Transmitter abhangt, von der Ausrichtung der Einzelneuronen 
gegenuber der Transmitterachse, und von weiteren physiologischen Parametern, 
sowie von der Art des gesendeten Signals. Ubersetzungsoperatoren sind 
Nichtobservable. 

20 Es ist bekannt, dass die genannten Nichtobservablen zwischen verschiedenen 
Personen urn Faktoren von bis zu 100 verschieden sind, und selbst innerhalb 
verschiedener Gebiete eines Gehims in vivo sehr verschieden sein konnen, und 
uberdies (im Vergieich zu EM Oszillationen langsam) schwanken konnen. Die 
resultierende Vielfalt theoretisch und praktisch moglicher Moden ist daher nahezu 

25 beliebig. Kontrollierte lokale Modifikation ohne lokale Kalibrierung ist daher im 
Normalfall nicht moglich, andererseits ist unter der Annahme relativ konstanten 
endogenen Inputs translokale Kalibrierung entbehrlich. 

S 800 (Spezifizieren der Mindest-Kalibrierungsvoraussetzungen): 

30 r , a, b, c, d, p x , p yt Ubersetzungsoperatoren fur ein Sensor-Transmitter-Paar. 

Es ist lediglich lokales Kalibrieren erforderlich. 

D 1100 (Mtndest-Kalibrierungsvoraussetzungen): 

Lokales Kalibrieren von r, a, b, c, d, p x , p yi Ubersetzungsoperatoren fur ein 
35 Sensor-Transmitter-Paar. 

E 100 (Kalibrierung erforderlich?): 

Lokales Kalibrieren ist erforderlich. 

40 D 41 00 (Apparaturbeschrankungen): 
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Hangt von der verwendeten Apparatur ab, beispielsweise auf der Eingangsseite 
gewohnliche EEG-Klebeelektroden, Digitalisierung mit einer Sampling-Rate von 
200/sek., auf der Ausgangsseite beispielsweise eine Spule mit maximalem 
Magnetfeld von 0.5 Tesla. Die Beschrankung auf einen Transmitter bringt es mit 
5 sich, dass verschiedene Moden nicht unabhangig voneinander beeinflusst 
werden konnen. 

0 4200 (Gesundheitliche Grenzwerte fur exogene elektromagnetische 
Felder): 

1 0 Beinhaltet beispielsweise die Empfehlungen von [14]. 
S 2000 (Kalibrieren): 

Wird beim Durchlauf durch die Flussdiagramme (Figuren 21-24) naher 
ausgefuhrt. 

15 

D 1200 (Spezifisches Hirnaktivitatsmodell): 

D 1000 mit empirisch fur jede relevanten Hirnaktivitat des jeweiligen Anwenders 
im betrefFenden Zeitrahmen ermittelten r , a, b, c, d, p Xtend0 > P y ,endo> some D 1300. 

20 D 1300 (Quantifizierung des Einflusses exogener EM Felder auf 
Hirnaktivitat): 

Empirisch fur jede relevante Hirnaktivitat des jeweiligen Anwenders im 
betreffenden Zeitrahmen ermittelter Einfluss exogener EM Felder auf die 
betreffende Hirnaktivitat, dargestellt durch Kontrollvariable p xexo und p yexo , in den 

25 vereinfachten Wilson-Cowan-Gleichungen zu dieser Hirnaktivitat. 

S 3000 (Steuern oder Regeln von Hirnaktivitat mit Hilfe gezielter 
Veranderung von Kontrollvariablen sowie Messung/ Berechnung von HAC): 

Wird beim Durchlauf durch das betreffende Flussdiagramm (Figur 25) genau 
30 angegeben. 

Figur 21: 

D 1100 (Mindest-Kalibrierungsanforderungen): 

35 Wie vorher. 

D 1110 (Kalibrierungsart): 

Lokal. 



40 D 1120 (Kalibrierungsumfang): 

Sensor S 1f Transmitter Ti. 
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E100 (Kalibrierung?): 

Ja. 

5 S 21 00 (Lokales Kalibrieren): 

Wird beim Durchlauf durch die Flussdiagramme (Figuren 22-23) naher 
ausgefuhrt. Im vorliegenden Beispiel gibt es nur ein einziges Sensor-Transmitter- 
Paar, fur das kalibriert werden muss. 

1 0 S 2300 (Identifizieren identischer Hlmaktivitaten): 

Entfallt, da verschiedene Sensoren und Transmitter in diesem Beispiel nicht zur 
Verfugung stehen. 

E 120 (Translokale Kalibrierung?): 

15 Nein. 

D 1200 (Spezifisches Hirnaktivitatsmodell): 

Wie vorher. 

20 D 1300 (Quantifizierung des Einfiusses exogener EM Feider auf 
Hirnaktivitat): 

Wie vorher, 

S 3000: 

25 Siehe Figur 25. 
Figur22: 

S 2110 (Messen EM Hirnaktivitat mit Si): 

Mit Hilfe der linksfrontalen Elektrode sowie der Referenzelektrode werden 
30 Potentialdifferenzen erfasst, und mit 100 Werten pro Sekunde digitalisiert. Die 
entstehende Zeitreihe (D 2205) wird gespeichert. 

D 2205 (Empirische Zeitreihe fur Si): 

Die aus S 21 10 entstehende Zeitreihe von Potentialdifferenzen. 

35 

S 2120 (Berechnen von Ist-Verlaufen von HAC aus der Zeitreihe von Si): 
Mit einer Sampling-Rate von 100/sek. erhalt man Frequenzen zwischen 1 und 50 
Hz. Per FFT (Fast Fourier Transformation) erhalt man fur ein Fenster 
festgelegter Lange (100 Messpunkte) das Leistungsspektrum. Hieraus erfolgt 
40 das Berechnen von Pos gemaft Formel in S 500. (Bereits bei HAC(100) beginnt 
hierS2130.) 
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D 2210 (1st -Verlaufe von HAC fur Si): 

Ist-Verlauf des Leistungsspektrums, also HAC(100), HAC(101), HAC(102), etc. 
Zusatzlich Ist-Verlauf von Pos, d.h. Pos(100), Pos(101), Pos(102), etc. 

5 

D 1000 (Generisches Hirnmodell): 

Siehe vorher. 

D 1049 (Einflussannahmen): 

10 Genauer erlautert bei D 1050, Figur 23. 

D 1051 (Datenbank): 

Die unendlich vielen im generischen Hirnmodell moglichen 
Parameterkombinationen werden a priori auf rechnerisch bewaltigbar viele 

1 5 reduziert: erstens durch Ausrechnen der Grenzen des Parameterbereichs, in dem 
endogene Inputs existieren, die zu einem nichtkonstanten EM Output fiihren, 
zweitens durch Diskretisieren eines n-fachen dieses Bereichs (Definieren eines 
Parametergitters, n naturliche Zahl). Jeder Parametersatz ist in diesem Gitter 
enthalten. Alle endogenen Inputs sind hier (einfachster Fall) als konstant 

20 betrachtet. Als Startwerte fur die Berechnung der Aktivitatsvariablen x und y in 
den vereinfachten Wilson-Cowan-Gleichungen werden im Intervall [0,1] 
gleichverteilte Zufallsvariable gewahlt. 

Es wird nun exemplarisch, ohne Anspruch auf numerische Prazision (bei vorerst 
auf Null gesetztem exogenen Input) vorgestellt, welchen Auswirkungen 

25 unterschiedliche Auspragungen der Parameter bei gleichem (konstantem) 
endogenem Input auf das Leistungsspektrum der resultierenden Mode haben 
(Fur die Figuren 7 bis 19 gilt: jeweilige x-Achse: Frequenz in Hertz, jeweilige y- 
Achse: quadrierter Absolutbetrag des zugehdrigen Fourierkoeffizienten): 
Figur 7 zeigt das HAC von Moden unter Variation der Membrankonstanten r : 

30 (i) wird erzeugt durch r = 5 , a=22, b=20, c=14, d=3, bei p x = 1.5 und p y = 0 . 

(ii) wird erzeugt durch r = 10, a=22, b=20, c=14, d=3, bei p x =1.5 und p y =0. 

(iii) wird erzeugt durch r = 15, a=22, b=20, c=14, d=3, bei p x =1.5 und p y =0. 

(iv) wird erzeugt durch r = 20 , a =22, b=20, c=14, d=3, bei p x = 1.5 und p y =0. 

Diese Moden sind unter konstantem Input stationar, sodass das gezeigte HAC, 
35 aufgetragen uber die Zeit, auch der HAC-Verlauf ist. Wenige Millisekunden 
dauernde Transiente, die von den Startwerten in die Oszillation fiihrten, wurden 
grundsatzlich weggelassen. 

Figur 8 zeigt das HAC von Moden unter Variation des 

Selbsterregungsparameters a: 
40 (i) wird erzeugt durch r = 10, a=14, b=20, c=14, d=3, bei p x =1.5 und p y =Q. 

(ii) wird erzeugt durch x=1 0, a=1 8, b=20, c=1 4, d=3, bei p x = 1 .5 und p y = 0. 
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(iii) wird erzeugt durch r = 10, a=22, b=20, c=14, d=3, bei p x =1.5 und p y =0. 

(iv) wird erzeugt durch r = 10 , a=26, b=20, c=14, d=3, bei p x =1.5 und /?, = 0 . 

Deutlich zu erkennen ist fehlende Oszillation bei den beiden kleineren 
Selbsterregungswerten. 
5 Figur 9 zeigt das HAC von Moden unter Variation des 
Fremdhemmungsparameters b: 

(i) wird erzeugt durch t = 10 , a =22, b=1 5, c=1 4, d=3, bei p x = 1 .5 und p y =0. 

(ii) wird erzeugt durch r = 10, a=22, b=20, c=14, d=3, bei p x =l.S und p y =0. 
(Hi) wird erzeugt durch r = 10, a=22, b=25, c=14, d=3, bei p x =1.5 und p y =0. 

1 0 (iv) wird erzeugt durch t = 10 , a=22, b=30, c=14, d=3, bei p x = 1.5 und p y =0. 

Figur 10 zeigt das HAC von Moden unter Variation des 
Fremderregungsparameters c: 

(i) wird erzeugt durch r = 10, a=22, b=20, c=10, d=3, bei p x =1.5 und p y =0. 

(ii) wird erzeugt durch r = 10, a =22, b=20, c=14, d=3, bei p x =1.5 und p y =0. 
15 (iii) wird erzeugt durch r = 10, a=22, b=20, c=18, d=3, bei p x =1.5 und p y = 0. 

(iv) wird erzeugt durch r = 10, a=22, b=20, c=22, d=3, bei p x =1.5 und p y =0. 

Deutlich zu erkennen ist fehlende Oszillation bei starker Fremderregung der 
inhibitorischen Neuronen innerhalb des Ensembles (Fall iv). 
Figur 11 zeigt das HAC von Moden unter Variation des 
20 Selbsthemmungsparameters d: 

(i) wird erzeugt durch r = 10, a =22, b=20, c=1 4, d=0, bei p x = 1 .5 und p y = 0 . 

(ii) wird erzeugt durch r =10, a =22, b=20, c=14, d=3, bei p x =1.5 und p y =0. 

(iii) wird erzeugt durch r = 10, a=22, b=20, c=14, d=6, bei p x =1.5 und p y =0. 

(iv) wird erzeugt durch r = 10, a=22, b=20, c=14, d=9, bei p x =1.5 und p y =0. 
25 Deutlich zu erkennen ist fehlende Oszillation bei starker Selbsthemmung (Fall iv). 

S 2115 (Fur jeden Satz von Nichtobservablen: Berechnen der zugehdrigen 
Modellmode und/oder Auffinden in der Datenbank): 

Einsetzen in die vereinfachten Wilson-Cowan-Gleichungen und numerisches 
30 Losen derselben. 

S 2120 (Berechnen von HAC-Verlaufen aus den Modellmoden und/oder 
Auffinden in der Datenbank): 

HAC-Berechnung auf die Ergebnisse von S 2115 anwenden. 

35 

D 1030 (HAC-Modellmoden): 

Ergebnisse von S 2120, auszugsweise dargestellt unter D 1051. 



D 1040 (HAC-Messinvarianten und -annahmen): 
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Um in empirischen Moden Modellmoden wledererkennen zu konnen, werden im 
einfachsten Fall HAC gewahlt, die ganz oder teilweise gegenuber der gewahrten 
Messmethode invariant sind, oder fur die diese Invarianz angenommen werden 
kann. Im vorliegenden Beispiel (Leistungsspektrum) ist diese Annahme: die 
5 Abschwachung zwischen dem Trager einer Modellmode und dem Sensor ist 
frequenzunabhangig. Das heiftt, das Verhaltnis der Koefflzienten einer 
Modellmode bleibt auf dem Weg zum Sensor erhalten. Als weitere Annahme wird 
das Signal als stuckweise stationar angesehen (mit Stationaritat im Mittel deutlich 
langer ais der Lange des Fourier-Fensters). Weitere Annahme es (s. [1]), dass 
10 die Potentialdifferenzen der an der betreffenden Hirnaktivitat beteiligten 
erregenden und die der hemmenden Neuronen zu gleichen Teilen in das EEG- 
Signal eingehen, und zwar als x - y (d.h. die Mischungswinkel sind a = 1 , p =-1). 

S 2130 (Fur 3i: Dekomposition von Ist-Verlaufen von modellnahen HAC in 
15 HAC-Modellmoden, mit Relevanzprufung): 

Grundsatzlich sind als modellnahe HAC die EM GroBen voreingestellt, die im 
Hirnaktivitatsmodell beschrieben werden (bspw. Potentialdifferenzen). Im 
vorliegenden Beispiel (quasistationare Neuraloszillatoren) wird stattdessen das 
Leistungsspektrum verwendet. Die a-priori-Anderung der Voreinstellung je nach 

20 verwendetem Hirnaktivitatsmodell ist dem Fachmann uberlassen. Hier Auffinden 
von Modellmoden, deren HAC-Modellmode, im einfachsten Fall dividiert durch 
einen Abschwachungsfaktor, bezuglich einer Norm (beispielsweise Li) 
bestmoglich das Leistungsspektrum erklart. Entfemen der HAC- 
Modellmode/Abschwachungsfaktor aus dem Leistungsspektrum. Auffinden 

25 weiterer Modellmoden, deren HAC-Modellmode bezuglich der Norm bestmoglich 
das neue Leistungsspektrum erklart, etc., bis ein vordefinierter Schwellenwert 
unterschritten wird. 

Wenn zum Beispiel das Ist-Leistungsspektrum modulo Rauschen stationar 
Maxima bei 10 Hertz mit Betrag 10, bei 20 Hertz mit Betrag 5, sowie kleinere 

30 Maxima bei weiteren Vielfachen von 10 Hz aufweist, so erklaren in erster 
Naherung (Maxima der Peaks statt Flache unter den Peaks) beide in Figur 12 
dargestellten HAC-Modellmoden das Ist-Leistungsspektaim annahemd gleich 
gut, das erste (E^) mit einer Abschwachung von ca. 10 / 2.2, das zweite (E 1f2 ) 
mit einer Abschwachung von ca. 10 / 6.8. Entfemen bspw. der ersten HAC- 

35 Modellmode*0.22 liefert annahemd das Nullspektrum. Analog fur die zweite HAC- 
Modellmode*0.68. Das heilit, im dargestellten Beispiel sind lediglich die HAC- 
Modellmoden Ei.i und Ei, 2 feststellbar. Zwischen diesen Moden kann man bei 
gleichbleibenden endogenen Input nicht unterscheiden. Daher erfolgt in diesem 
Falle die Unterscheidung mit Hilfe von Testsignalen. 

40 
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D 2300 (Liste relevanter HAC-Modellmoden im Messbereich jedes Sensors 

Si): 

Besteht hier nur aus E if1 und Ei, 2 , beide fur den Sensor Si. 

5 D 2350 (Liste moglicher Modellmoden zu jedem Element von D 2300): 

Diese Liste enthalt, neben anderen, die Modellmoden Mi.i.i (fur Ei.O und Mi i2 ,i 
(fur Ei, 2 ), deren Erzeugung im folgenden behandelt wind. 

D 2400 (Liste moglicher Parametersatze mit endogenem Input zu jeder 

10 Modellmode in D 2350): 

Diese Liste enthalt, neben anderen, zu jeder Modellmode aus D 2350 
Modellmoden Parametersatze und endogene Inputs, die diese Modellmode 
gemaB der vereinfachten Wilson-Cowan-Gleichungen erzeugen. D 2400 umfasst 
hier bspw. fur Mi l1f1 den Parametersatz r = 10, a=22, b=20, c=14, d=3, bei 

15 endogenem Input p x = 1.5 und p y = 0, sowie fur Mi, 2 ,i den Parametersatz r =10, 

a=26, b=20, c=18, d=3 bei endogenem Input ^=1 .36 und ^=-0.14. 
Fiaur 23: 

20 D 2300, D 2350, D 2400, D 1000, D 1051, D 4100, D 4200: 

Wie vorher. 

D 1050 (Einflussinvarianten und Annahmen): 

Die Kontrollvariablen in den Gleichungen des generischen Hirnmodells entstehen 
25 durch Anwendung des zu einer Modellmode und zum verwendeten Transmitter 
^ gehdrigen (noch unbekannten) Ubersetzungsoperators 0 auf ein 
Transmittersignal. Im vorliegenden Beispiel: 

rx = -x + S(ax -by + p xendo + U x (Transmittersignal)) 

ry = -y + S(cx -dy + p ymdo + U y (Transmittersignal)) 

30 Bezuglich der Gestalt des Ubersetzungsoperators werden (uberpriifbare) 
Annahmen gemacht, z.B. fur einen Spezialfall akustischer Stimulation 0 X = Cl y (s. 
[18]). Die gleiche Annahme wlrd hier getroffen, zusatzlich die Annahme, dass im 
Mittel, innerhalb eines Geltungsbereichs, zwei sich lediglich bezuglich ihrer 
Amplitude unterscheidende Transmittersignale (bspw. Amplitudenverhaltnis 1:r, r 

35 reelle Zahl >0) sich in der Amplitude ihrer zugehorigen exogenen Inputs ebenfalls 
wie 1:r verhalten. Weitere Annahmen sind: p xemh , p y mh wahrend der lokalen 

Eichung stationar. 
S 2160, D 2500: 
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Als Testsignal kommt zunachst jedes beliebige, im Rahmen der 
Apparaturbeschrankungen (D 4100) und Gesundheitsgrenzwerte (D 4200) 
realisierbare Signal in Frage. 

Aufgabe des Testsignals ist es erstens, aus D 2400 Kandidaten fur erzeugende 
5 Parameter und endogene Inputs der Modellmoden der beobachteten HAC- 
Modellmoden solange zu eliminieren, bis deren Zuordnung zur jeweiligen 
beobachteten HAC-Modellmode eindeutig ist, und zweitens die Ubersetzung des 
Testsignals in exogenen Input festzustellen. 

Aus der Menge aller moglichen Testsignale konnen Testsignale, die diese 
10 Bedingungen erfullen, numerisch ermittelt werden (bspw. Wahl von geordneten 
Funktionenbasen plus brute force - Berechnung). 

Eine zulassige Vereinfachung ist es, von einem Testsignal sowohl einen Bezug 
zu den gemessenen Hirnaktivitaten zu verlangen, als auch zum verwendeten 
Hirnaktivitatsmodell, als auch in selbstkonsistenter Weise zu einem bei der 

15 Modifikation in S 3000 einzusetzenden Signal. Bei linearen Oszillationen wird 
zwischen zur Oszillationsfrequenz inkommensurablen sowie kommensurablen 
Frequenzen unterschieden. Im vorliegenden Beispiel wird analog fur die 
Grundfrequenz der Oszillation (kleinste im Leistungsspektrum enthaltene, von 
den Frequenzen eines Zufallssignals unterscheidbare Frequenz) vorgegangen: 

20 Nach Figur 12 ist die Grundfrequenz der beiden ahnlichen HAC-Modellmoden 10 
Hz. Als inkommensurable Testsignale kommen somit sinusformige Signale der 
Frequenzen 3 Hz, 7 Hz, 9 Hz, 11 Hz etc. in Frage. Als kommensurable 
Testsignale kommen sinusformige Signale der Frequenzen 2 Hz, 4 Hz, 5 Hz, 6 
Hz, 8 Hz, 10 Hz etc. in Frage. (Phasenverschiebungen sind im vorliegenden 

25 Beispiel fur das Leistungsspektrum einer Mode auf einem Trager irrelevant.) Des 
weiteren sind grundsatzlich Elementarloops fur Testsignale verwendbar, d.h. 
zeitverzogerte Ruckkopplungen eines Messsignals, beispielsweise der Gestalt 
exo(t)=x(t-10)-y(t-10). 

Neben die dargestellten Moglichkeiten der Gestalt des Testsignals treten 
30 verschiedene Amplitudenmultiplikatoren, beginnend mit einer geringen Amplitude. 
Beispiel: Testsignal =(0-05 Tesla, sin(2*Pi*t*3/1000), (1,2,3,4)) , t Zeit in 
Millisekunden. 

S 2170 (Senden des Testsignals): 

35 Erfolgt in der in D 2500 festgelegten Reihenfolge, mit Pausen zwischen je zwei 
Sendevorgangen in der Grofcenordnung von Vielfachen der 
Membranzeitkonstante, zur Vermeidung temporaler Summationseffekte. 

S 2110, D 2520: 

40 Wie vorher. Falls Messung und Senden des Testsignals gleichzeitig erfolgen, und 
beispielsweise fehlende oder unzureichende Sensorabschirmung zu Artefakten 
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fiihrt, so sind diese aus dem Signal herauszurechnen. Die Berechnung der 
Starke des von einem Transmitter in einem Sensor erzeugten Artefakts erfolgt 
ublicherweise uber die Entfernung zwischen beiden, sowie Materialgleichungen 
der zwischen beiden befindlichen Substanzen. Zusatzlich besteht die Moglichkeit, 
5 diese Artefakte vor Beginn einer Kalibrierung direkt zu messen, indem von 
naturlichen Aktivitaten gut abgrenzbare Signale sehr geringer Amplitude 
gesendet werden (so gering, dass ein Einfluss auf neuronale Aktivitaten als 
unwahrscheinlich angenommen werden darf). Ohne Transmittersignal moglicher 
Artefakte (bspw. durch Muskelbewegungen wahrend der Messung) werden 
10 konventionell beseitigt 

S2120: 

Wie vorher. 

15 S 2135 (Dekomposition nach relevanten HAC-Modellmoden mit Testsignal): 

Im Gegensatz zu S 2130 wird hier nicht nach beliebigen HAC-Modellmoden 
zerlegt, sondern nach HAC-Modellmoden mit Testsignal. 
Hierbei sind lediglich Parametersatze und endogener Input aus D 2400 
zugelassen, die jetzt durch exogenen Input zu ggf. anderen Losungen 
20 (Modellmoden mit Testsignal), und damit zu ggf. anderen HAC-Modellmoden mit 
Testsignal fiihren. 

Die Abschwachungen von S 2130 bleiben erhalten, wenn angenommen wird, 
dass sich an der Lokalisierung des Tragers der jeweiligen Mode durch exogenen 
Input nichts andert. 

25 Die zu berucksichtigenden HAC-Modellmoden haben beispielsweise unter 
sinusformigem exogenem Input der Frequenz 7 Hertz fur Multiplikatoren 2,4,6,8 
die fur Mode I (siehe Figur 12 (i)) in Figur 13 dargestellte Gestalt: 
(0 fiir Pexo =2*sin(2*tf *7/1000) 
(ii)ftir p m = 4*sin(2/rt*7/1000) 

30 (iii) fur p ao =6*sin(2^*7/1000) 
(iv) fur p m = 8*sin(2***7/1000) 

Hieraus ergeben sich bereits zur Veranschaulichung der Quantifizierung des 
Einflusses des Signals einige Beobachtungen: Das bezuglich der 7 Hertz und der 
14 Hertz Komponente invertierte Leistungsspektrum normalisiert sich beim 
35 Verdoppeln des Multiplikators von 2 auf 4. Der quadrierte Absolutbetrag des 7 
Hertz - Anteils weist bei der Erhohung des Multiplikators von 2 auf 4 folgendes 
(ungefahre) Verhaltnis auf: 1/4, von 4 auf 6: 3/2, von 6 auf 8: 4/3. 

Der selbe exogene Input erzeugt im Falle der anderen HAC-Modellmode (Mode 
40 II, siehe (ii) aus Figur 12) die in Figur 14 fur Multiplikatoren 2,4,6,8 dargestellten 
Leistungsspektren: 



WO 03/085546 



-33- 



PCT/EP03/03545 



(i) fur p m = 2*sin(2;r**7/1000) 

(ii) fCir Pexo =4*sin(2^*7/1000) 

(iii) fur p m =6*sin(2^*7/1000) 

(iv) fur p m =8*sin(2^*7/1000) 

5 Bereits bei kleinen Multiplikatoren trennt das Testsignal deutlich zwischen Mode I 
und Mode II, wie aus dem Vergleich von Fig. 13 (i) und Fig. 14 (i) ersichtlich ist. 

Das Vervielfachen der Signalintensitat h gestattet die Einflusskalibrierung. Wenn 
beispielsweise im Fall von Mode II das Verdoppeln der Signalintensitat h auf l 2 
10 das Verhaltnis von erstem Peak zu zweitem Peak von 6 auf 2.5 absinken lasst, 
so ergibt der Abgleich mit der HAC-Modellmodendatenbank beispielsweise, dass 
h einem Multiplikator von 2 entspricht (Fig. 14 (i) und l 2 einem Multiplikator von 4 
(Fig. 14(H)). 

15 Das vorangegangene Beispiel vermied aus Griinden der Darstellungseinfachheit 
die in der Realitat haufig auftretende Komplikation, dass Messen und Senden 
nicht gleichzeitig erfolgen konnen. Bei im Vergleich zur Oszillationsperiode 
kleinen Messunterbrechungen wird im Normalfall das ununterbrochene 
Messsignal mit ublichen signaltheorethischen Methoden rekonstruiert. Im 

20 folgenden wird stets en solches ununterbrochenes Messsignal als gegeben 
betrachtet. Transmissionsunterbrechungen werden beispielhaft fur den Fall von 
Elementarloops aus exogenem Input dargestellt. 

Als Beispiel fur ein Elementarloop wird das Messsignal x(t)-y(t) zeitverzogert mit 
25 einer Verzogerung T ruckgekoppelt, das heiftt, es ergeben sich beispielsweise 
fur eine Verzogerung von 10 Millisekunden die Gleichungen 
Tx = -x + S(ax-by + p xendo +p Xiexo ) und 

ry = ~y + S(cx -dy+ p yendo + p y exo ) >wobei 

Px y exo = Py,exo = M* ( x (t ""10)- y (t — 1 0)) ( we nn man davon ausgeht, 

30 dass der Ubersetzungsoperator das Signal nicht verzerrt und/oder ebenfalls 
verzogert). M ist hierbei ein gegebenenfalls uber sukzessive 

Amplitudenverdopplung des Transmittersignals zu ermittelnder Multiplikator, alle 
anderen Symbole wurden bereits erlautert. 

Urn den Fall nicht gleichzeitig aktiver Sensoren und Transmitter (relevant 
35 beispielsweise bei Blockierung des Sensors wahrend des Sendens) zu 

behandeln, wahlt man als Umschaltzeit beispielsweise die Verzogerungszeit. Im 
vorliegenden Beispiel schaltet man also alle 10 Millisekunden zwischen Messen 
und Senden hin -und her. In den obigen Gleichungen ist dann: 



40 p Xl exo = Py.exo = H • Xi(t) mit 
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xi(t):=tf(t mod20 > 10, x(t-10)-y(t-10), 0) in generischem Code 

Figur 15 zeigt Mode I mit diesem exogenen Input fur ju = 4 (i) und ju = 8 (ii). 
5 Figur 16 zeigt Mode II mit diesem exogenen Input fur ju = 4 (i) und // = 8 (ii). 
Figur 17 zeigt Mode I mit diesem exogenen Input fur // = 12 (Hi) und // = 16 (iv). 
Figur 18 zeigt Mode II mit diesem exogenen Input fur n = 12 (iii) und // = 16 (iv). 

Auch hiertrennt das Testsignal zwischen den Moden I und II, beispielsweise iiber 
10 Frequenzverschiebung (Figur 15 (ii) gegenuber Figur 16 (ii)) und/oder dem 

Verhaltnis der Maxima aufeinanderfolgender Peaks (Figur 17 (iii) gegenuber 

Figur 18 (iii)). Das Verhaltnis der Maxima aufeinanderfolgender Peaks bei 

verschiedenen exogenen Inputs in eine Mode gestattet auch hier, analog zum 

dargestellten Sinus-Input, die Einflusskalibrierung. 
15 Phasenbetrachtungen unterbleiben in diesem einfachen Beispiel. Bei einer Fine- 

Tuning-Kalibrierung .werden zusatzlich Sequenzen phasenverschobener 

Testsignale verwendet. 

Fur die weitere Darstellung von Beispiel 1 wird davon ausgegangen, dass mit 
20 Hilfe zeitverzogert riickgekoppelter Testsignale festgestellt wurde, dass der 
beobachteten Hirnaktivitat die Mode I unterliegt, mit ihren Nichtobservablen 
r = 10, a=22, b=20, c=14, d=3, bei p Xitndo =L5 und />, >CTrfo =0, in die das 

dargestellte zeitverzogert ruckgekoppelte Hertz Signal des Transmitters Ti mit 
40*[Transmitterfeldstarke in Tesla]*xi(t) eingeht, und die auf dem Weg von ihrem 
25 Trager zum Sensor Si einer Abschwachung von 1 0/2.2 unterliegt. 

Anmerkung 1:' Im allgemeinen ist nicht unbedingt davon auszugehen, dass 
geeignete Moden durch Parameter und Inputs des gewahlten Gitters erzeugt 
werden, d.h. es sind lokale Verfeinerungen des Gitters und/oder Interpolationen 
30 durchzufuhren. 

Anmerkung 2: Inkonsistente Ergebnisse innerhalb einer Signalsequenz weisen 
darauf hin, dass wahrend des Sendens ein oder mehrere Sprunge des 
quasistationaren endogenen Inputs stattgefunden haben. In diesem Fall mussen 
Moden mit zeitveranderlichen Input fur die Kalibrierung berucksichtigt werden. 

35 

E 210 (Mehrere bestapproximierende Satze von Parametern): 

In diesem Beispiel nicht. Eindeutigkeit ist bis auf Aquivalenz definiert. Zwei Satze 
A und B von Parametern heilXen aquivalent, wenn unter alien Inputs die jeweilige 
A-HAC-Mode von der B-HAC-Mode innerhalb einer vordefinierten Fehlerschranke 
40 ununterscheidbar ist. 
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Fiaur24: 

Entfallt, da in diesem Beispiel nur lokale Kalibrierung benotigt wird. 
Fiaur25: 

5 

D 2100, D 1200, D 1300, D 1050, 0 1051: 

Wie vorher. 

D 1010 (HAC-lst-Verlaufe): 

10 Aus den laufenden Messungen von Himaktivitat wird Pos laufend berechnet 

30 12 

CLaufend" := „ in geeigneten Zeitabstanden"). pos := ^fi l^fi gemafc S 500. 

i=I3 /=8 

S 2110, S 2120: 

Wie vorher. 

15 

D 4100, D 4200: 

Wie vorher 

D4300: 

20 Kann mit nur einem Transmitter nicht realisiert werden. Wird in Beispiel 3 
erlautert. 

S 3100 (Berechnen eines oder mehrer Transmittersignale): 

Die in der Regel unendliche Vielzahl moglicher Transmittersignale wird analog 
25 zur Berechnung der Testsignale reduziert. Im einfachsten Fall werden 
Transmittersignale verwendet, die bereits als Testsignale verwendet wurden. Die 
Transmittersignale werden auf Machbarkeit (D 4200, D 4100, D 4300) gepriift 
und bezuglich vordefinierter Nutzenfunktionen geordnet (beispw. moglichst 
niederfrequent, moglichst kleine Feldstarke, etc.). Das bestmdgliche 
30 Transmittersignal wird dann automatisch selektiert. 

Beispiel: 0.4 Tesla oben genanntes zeitverzogert riickgekoppeltes Signal des 
Transmitters Tv fiihrt zum Erreichen des Verhaltenziels, da Pos mehr als 
verdoppelt wird (vgl. Fig. 12 (i) mit Fig. 17 (iv)). 

35 S 3400 (Senden): 

Senden des/der gemali S 3100 ausgewahlten Transmittersignale (hien mit dem 
Transmitter Ti). 

D 4040 (HAC-Verlaufsprognose): 
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Die Modellmoden mit exogenen Inputs bestirnmen die entsprechende 
Hirnaktivitat im Zeitverlauf, und liefern eine Zeitreihe, aus deren in der Zukunft 
liegenden Teilen wie vorher eine theoretische HAC-Zeitreihe berechnet wird, 
letztere ist die HAC-Veriaufsprognose. Im vorliegenden Beispiel ist nur eine 
5 Modellmode relevant, die HAC-Verlaufsprognose ergibt sich daher durch Teilen 
des in Figur 13 (I) dargestellten HAC durch den Abschwachungsfaktor 10/2.2. 

S 3500 (Vergleichen): 

Zum einen gilt hier: Subtraktion des Ist-HAC-Verlaufs aus S 2120 und des 
10 prognostizierten HAC-Verlaufs (D 4040), Betragsbildung fur jede Frequenz 
zwischen 1 und 50 Hertz (Abweichung = Maximum dieser Betrage), zum anderen 
wird das Erreichen des Verhaltensziels (des Soll-HAC-Verlaufs) analog durch 
Vergleichen von Ist- und Soll-HAC-Verlauf ermittelt. 

15 E 300 (Abweichung akzeptabel?): 

Treffen die Annahmen und die Ergebnisse der Kalibrierung zu, so ist davon 
auszugehen, dass sich Prognose- und Ist-HAC nur unwesentlich unterscheiden, 
(Abweichung < vordefinierter Schwellenwert) man kann demzufolge mit dem 
Senden fortfahren (S 3400). Ist dieses nicht der Fall, so gilt es per numerischer 
20 Losung der Gleichungen des Hirnmodells mit anderen endogenen Inputs (S 
3600) zu prtifen, inwieweit Sprunge des stuckweise stationaren endogenen 
Inputs die Abweichung erklaren. 

Im Falle des Nichterreichens des Verhaltensziels innerhalb der berechneten 
25 Zeitspanne erfolgt der Ubergang zu S 3800, anderenfalls die Meldung 
„Verhaltensziel erreicht". 

S 2000 (Kalibrieren) 

Wie bereits erlautert. 

30 

Abschliefcend ist zu diesem Beispiel zu bemerken, dass das Hirnaktivitatsziel 
durch Kalibrieren und anschlieftende kontrollierte Modifikation zuverlassig 
erreicht wird, und somit bei Zutreffen des reduzierten Davidson-Modells auch das 
Verhaltensziel, 

35 

Beispiel 2 (Schritte werden nur auszuqsweise erlautert): 



40 



Fiaur 20: 

D 2000 (Verhaltensziel): 
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Verstarken positiver Emotionen und Abschwachen negativer Emotionen fur die 
Dauer der Anwendung 

D 3000 (Verhaltensmodell): 

5 Davidson-Modell (siehe z.B. [1 1]), wonach positive Emotionen mit einem hoheren 
Quotienten von Beta (13-30 Hz) - EEG-Aktivitat zu Alpha (8-12 Hz) - EEG- 
Aktivitat im linksfrontalen Kortex einhergehen, und negative Emotionen mit einem 
hoheren Quotienten von Beta (13-30 Hz) - EEG-Aktivitat zu Alpha (8-12 Hz) - 
EEG-Aktivitat im rechtsfrontalen Kortex. 

10 

S 500 (Spezifizieren HAC. Spezifizleren Soll-Verlaufe von HAC): 

30 12 

Analog zu Beispiel 1. Abgeleitete HAC sind Pos:=^f^/^f^ fur einen 

/=I3 ' /«8 

30 12 

linksfrontalen Sensor Si und Neg:=Y< f^'Y^f* fur einen rechtsfrontalen Sensor 

S 2 . (fy Koeffizient der i-Hertz-Mode des Sensors j). Pos soil erhoht werden, also 
15 beispielsweise Pos(mit Einfluss) > 2 *Pos(ohne Einfluss), Neg soli verringert 
werden, also beispielsweise Neg(mit Einfluss) < 0.5*Neg(ohne Einfluss). Diese 
Anforderung definieren die Soll-HAC Verlaufe. Berechnung wie in Beispiel 1 . 

D 4000 (Minimale Apparaturvoraussetzungen): 

20 Oberflachen-EEG mit einer Elektrode S<[ linksfrontal mit Referenzelektrode am 
linken Ohr, einer Elektrode S 2 rechtsfrontal mit Referenzelektrode am rechten 
Ohr (s. [16]), zwei extrakranialen Spulen (Transmitter Ti und T 2 ), wobei Ti 
beispielsweise auf einer an der Kopfoberflache verlaufenden Verbindungslinie 
zwischen Si und einer zentralen Elektrode Cz (s. [16]) in unmittelbarer Nahe von 

25 Si lokalisiert ist, und T 2 analog zwischen S 2 und Cz, wobei die Entfernung 
zwischen S^ und Ti gleich der Entfernung zwischen S 2 und T 2 ist. Alles andere 
ist wie in Beispiel 1 . 

D 2100 (Soll-Verlaufe von HAC): 

30 Pos (mit Einfluss) > 2*Pos (ohne Einfluss) und Neg (mit Einfluss) < 0.5*Neg 
(ohne Einfluss) 

D 2200 (HAC-Berechnungsvorschriften): 

Wie in Beispiel 1, zusatzlich fur den zweiten, rechtsfrontalen Sensor 

35 Neg :=£^ 2 /i^ 2 . 

D 1000 (Generisches Hirnaktivitatsmodell): 

Vereinfachtes Wilson-Cowan-Modell, wie in Beispiel 1. 
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S 800 (Spezifizieren der Mindest-Kalibrierungsvoraussetzungen): 

r, a, b, c, d, p x ,p y1 Ubersetzungsoperatoren fur die vier Sensor-Transmitter- 

Paare Si-Ti, Si-T 2 , S 2 -Ti, S 2 -T 2 . Wieder geniigt lokales Kalibrieren, da beide 
5 frontalen Hirnbereiche vereinfacht als unabhangig angenornmen werden, sowie 
der jeweilige endogene Input als konstant. 

D 1100 (Mindest-Kalibrierungsvoraussetzungen): 

Lokales Kalibrieren von t , a, b, c, d, p x ,p y , die Ubersetzungsoperatoren fiir die 
1 0 vier Sensor-Transmitter-Paare S r Ti , Si-T 2 , S 2 -Ti , S 2 -T 2 . 

E 100 (Kalibrierung erforderlich?): 

Lokales Kalibrieren ist erforderlich. 

15 D 41 00 (Apparaturbeschrankungen): 

Analog zu Beispiel 1, auf der Ausgangsseite beispielsweise zwei Spulen mit 
jeweils maximalem Magnetfeld von 0.5 Tesla. 

D 4200 (Gesundheitliche Grenzwerte fur exogene EM Felder): 

20 Wie in Beispiel 1. 

S 2000 (Kalibrieren): 

Siehe spater. 

25 D 1 200 (Spezifisches Hirnaktivitatsmodell): 

Wie in Beispiel 1 . 

D 1300 (Quantifizierung des Einflusses exogener EM Felder auf 
Hirnaktivitat): 

30 Wie in Beispiel 1 . 

S 3000 (Steuern oder Regeln von Hirnaktivitat mit Hilfe gezielter 
Veranderung von Kontrollvariablen sowie Messung/ Berechnung von HAC): 

Wird beim Durchlauf durch das betreffende Flussdiagramm (Figur 25) genau 
35 angegeben. 

Fiaur21: 

D 1100 (Mindest-Kalibrierungsanforderungen): 

40 Wie vorher. 
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D 1110 (Kalibrierungsart): 
Lokal. 

D 1120 (Kalibrierungsumfang): 

5 Sensoren Si , S 2i Transmitter Ti , T 2 . 

E 110 (Keine Kalibrierung?): 

Ja. 

10 S 2100 (Lokales Kalibrieren): 

Analog zu Beispiel 1, fur die vier Sensor-Transmitter-Paare S r Ti, S r T 2 , S^T^ 
S 2 -Ti. Wesentliche Vereinfachung beim vierfachen Durchlaufen der lokalen 
Kalibrierung ist es, dass Parameter und endogene Inputs von durch Si 
feststellbaren Moden bereits beim Durchlaufen der lokalen Kalibrierung fur das 

1 5 Paar S1-T1 festgestellt werden, und somit bei der lokalen Kalibrierung fur S r T 2 als 
gegeben angenommen werden (u.a. bei der Dekomposition). Analog fur Moden 
im Messbereich von S 2 . Es wird empfohlen, die lokale Kalibrierung eines Sensors 
mit demjenigen Transmitter zu beginnen, der gegenuber diesem Sensor eine 
minimale raumliche Distanz aufweist. 

20 Als Ergebnis der lokalen Kalibrierung ergibt sich beispielsweise: 
Fur Sensor Si: 

Mode I aus Beispiel 1 mit Einfluss von Ti als 40*[Transmitterfeldstarke von T1 in 
Tesla]*xii kaiibriert und Einfluss vn T 2 = 0. (xii(t) ist die zeitverzogerte 
Ruckkopplung des von Si gemessenen Signals, in rekonstruierter Form). 
25 Fur Sensor S 2 : 

Mode II aus Beispiel 1 mit Einfluss von T 2 als 40*[Transmitterfeldstarke von T 2 in 
Tesla]*xi 2 (t) kaiibriert und Einfluss von T1 = 0. (xi 2 (t) ist die zeitverzogerte 
Ruckkopplung des von S 2 gemessenen Signals, in rekonstruierter Form). 

30 S 2300 (Identische Hirnaktivitaten identifizieren): 

Obige Moden sind nicht identisch. 

E 120 (Translokales Kalibrieren?): 

Nein. 

35 

D 1200 (Spezifisches Hirnaktivitatsmodell): 

Wie vorher. 



D 1300 (Quantifizieren des Einflusses exogener elektromagnetischer Felder 
40 auf Hirnaktivitat): 

Wie vorher. 
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Fiauren 22 -23: 

Wie vorher, bzw. wie gerade unterS 2100 dargestellt. 

5 Fiaur 24: 

Enfallt, da nur lokale Kalibrierung benotigt wind. 

Fiaur 25: 

Analog zu Beispiel 1, bis auf S 3100, was im folgenden naher erlautert werden 
10 soli: 

Beispielsweise fuhrt das kombinierte Signal von Ti, xii(t) mit 0.4 Tesla und T 2 , 
xi 2 (t) mit 0.4 Tesla bei Mode I zu einer VergrdRerung von Pos auf uber 200% des 
ursprunglichen Wertes (siehe Beispiel 1), bei Mode II zu einer Verringerung von 
Neg auf unter 50% des ursprunglichen Wertes (vergleiche Figur 12 (ii) mit Figur 
15 18 (iv)). Das Hirnaktivitatsziel wird somit erreicht, und ebenso bei Zutreffen des 
Davidson-Modells auch das Verhaltensziel. 

Beispiel 3: 

20 Beispiel 3 soli nicht in der gleichen Ausfuhrlichkeit wie die vorangegangenen 
Beispiele (z.B. samtliche Moden, lokale Kalibrierung aller Sensor-Transmitter- 
Paare, etc) erlautert werden, sondem hauptsachlich beziiglich wesentlicher 
Abweichungen von den vorangegangenen Beispielen: 

25 Figur 20: 

D 2000 (Verhaltensziel): 

Verlangem mentaler Aufmerksamkeit fur die Dauer der Anwendung. 

30 D 3000 (Verhaltensmodell): 

Basierend bspw. auf [20] und [21], geht FmB (Frontal midline theta), also 6-7 HZ 
Aktivitat im Bereich der Fz-Elektrode (s. [16]) mit der Aufrechterhaltung 
konzentrierter Aufmerksamkeit bei mentaler Prozessen einher. Diese 
Feststellung wird im folgenden als Ishihara-Yoshii-Modell bezeichnet. 
35 Oberflachennahe Neuronen im Messbereich des Fz-Sensors sind mit einigen 
anderen Bereichen direkt verbunden, insbesondere wurde subkortikales Driving 
festgestellt. Leistungsspektrum und dessen Berechnung wie vorher, statt Pos 

Oder Neg diesmal fur jeden Sensor Att := Y^tf l^fi . (max=50) 

/=6 (=1 



40 D 4000 (Minimale Apparaturvoraussetzungen): 
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Zur Feststellung etwaiger Ausbereitungsphanomene ist es sinnvoll, die Fz- 
Elektrode durch zumindest eine benachbarte Elektrode mit physiologisch direkt 
mit dem Messbereich der Fz-Elektrode verbundenem Messbereich zu erganzen 
(beispielsweise der Cz-Elektrode), d.h. Oberflachen-EEG mit zwei Elektroden und 
5 zwei Referenzelektroden (wie in [21] dargestellt). Zwei Spulen als Transmitter 
sind jeweils in unmittelbarer Nahe der betreffenden Elektrode extrakranial 
lokalisiert: Ti anterior bezuglich Fz, T 2 posterior bezuglich Cz. Alles andere wie in 
Beispiel 1. 

10 D 21 00 (Soll-Verlaufe von HAC): 

Beispielsweise: Atti(mit Einfluss)= Atti(ohne Einfiuss, bei vollem Driving), Att 2 (mit 
Einfluss)<= Att 2 (ohne Einfiuss). Hierbei steht „1" fur „Fz und Referenzelektrode", 
„2" fur Cz und Referenzelektrode", , = M steht fur „die rechte Seite weicht nicht 
mehr als 20% von der linken Seite ab". 



D 2200 (HAC-Berechnungsvorschriften): 

Leistungsspektrum wie in Beispiel 1, zusatzlich fur jeden der zwei Sensoren 



20 D 1000 (Generisches Hirnaktivitatsmodell): 

Vereinfachtes Wilson-Cowan-Modell, wie in Beispiel 1, wobei jetzt die Rolle 
translokaler Effekte genauer betrachtet wird und die Konstanzannahme fur den 
endogenen Input aufgegeben wird. Der endogene Input in die Gesamtheit der 
erregenden Einzelneuronen des j-ten Oszillators lasst sich im einfachsten Fall der 



25 translokalen Kalibrierung aufspalten in r Xj ,endo\ / 3 J ■ 

Hierbei stellt w hj" die Starke des Einflusses des l-ten auf den j-ten neuralen 
Oszillator dar, und das betreffende A die Verzogerung des Wirksamwerdens 
dieses Einflusses. 

30 S 800 (Spezifizieren der Mindest-Kalibrierungsvoraussetzungen): 

Wie in Beispiel 2, plus Funktionalmatrix. Lokales und translokales Kalibrieren 
erforderlich. 

D 1100 (Mindest-Kalibrierungsvoraussetzungen): 

35 Wie in Beispiel 2, plus Funktionalmatrix. 



15 



/=6 /=1 




E 100 (Kalibrierung erforderlich?): 

Ja. 
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D 4100 (Apparaturbeschrankungen): 

Analog Beispiel 2. 

D 4200 (Gesundheitliche Grenzwerte fur exogene elektromagnetische 
5 Felder): 

Wie in Beispiel 2. 

S 2000 (Kalibrieren): 

Wird beim Durchlauf durch die Flussdiagramme (Figuren 21-24) naher 
10 ausgefuhrt. Es wild davon ausgegangen, dass die Stationaritat des endogenen 
Inputs zwar iiber die Kalibrierungsdauer, nicht aber iiber die Anwendungsdauer 
gegeben ist. 

D 1200 (Spezifisches Hirnaktivitatsmodell): 

15 Wie in Beispiel 1 . 

D 1300 (Quantifizieren des Einflusses exogener EM Felder auf Hirnaktivitat): 

Wie in Beispiel 1. 

20 S 3000 (Steuern Oder Regeln von Hirnaktivitat mit Hilfe gezielter 
Veranderung von Kontrollvariablen sowie Messung/ Berechnung von HAG): 

Wird beim Durchlauf durch das Flussdiagramm (Figur 25) genau angegeben. 

Figur 21 : 

25 

D 1100 (Mindest-Kalibrierungsanforderungen): 

Wie vorher. 

D 1 1 1 0 (Kalibrierungsart): 

30 Translokal (beinhaltet lokal). 

D 1120 (Kalibrierungsumfang): 

Wie in Beispiel 2. 

35 S 2100 (Lokales Kalibrieren): 

Wird beim Durchlauf durch die Flussdiagramme (Figuren 22-23) naher 
ausgefuhrt. 

S 2300 (Identische Hirnaktivitaten identifizieren): 

40 Analog zu Beispiel 2. 
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E 120 (Translokale Kalibrierung?): 

Ja 

S 2400 (Translokales Kalibrieren): 

5 Wird bei Figur 24 genau dargestellt. 

D 1200, D 1300: 

Wie in Beispiel 2, plus Funktionalmatrix. 
10 Figur 22: 

Analog zu Beispiel 2. Von Interesse ist jetzt bspw. im Messbereich von Sensor 
Si eine HAC-Modellmode mit unteriiegenden Parametern r=10, a=16, b=20, 
c=8, d=3, bei endogenem Input/?, = 1.5+m*sin(2**:*f*7/1000) und 

15 ^=1.5+w*sin(2*^*f*7/1000), dargestellt fur m=0 in Figur 19 (i), fur m=2 in 
Figur 19 (ii), und fur m=4 in Figur 19 (iii). (In diesem Beispiel ist angenommen, 
dass sich das subcortikale Driving mit 7 Hertz auf den erregenden und den 
hemmenden Teil des beschriebenen Neuronenensembles gleichermalien 
auswirkt.) Diese Mode soli als Mode III bezeichnet werden. Testsignale, die so 

20 nahe wie mdglich am 7 Hz sinusoidalen Driving sind, erscheinen als geeignet, 
beispielsweise ein 7 Hz Sinussignal, das alle 10 Millisekunden fur 5 Millisekunden 
abgeschaltet wird, urn die Messung zu gestatten. Ein solches Testsignal, mit 0.5 
Tesla vom Transmitter Ti abgestrahlt, soli beispw. mit dem Multiplikator 16 
wirksam werden, analog vom Transmitter T 2 abgestrahlt ohne Wirkung auf Mode 

25 III bleiben. 

Lokale Kalibrierung umfasst hier sinnvollerweise Phasenverschiebung von 
Testsignalen. Das HAC-Minimum der getriebenen Mode liefert die zum 
Ausschalten der Mode benotigte Phasenverschiebung 3 . Beispielsweise wird ein 
30 wie oben dargestellt unterbrochenes Sinussignal, das urn 3 verzogert mit 0.5 
Tesla vom Transmitter Ti abgestrahlt wird, die Mode III weitgehend 
ausschalten.(Ergebnis hiervon siehe Figur 19 i). 

Einzige S 2 -Mode ist beispielsweise die S r Mode aus Beispiel 1. Beide o.g. T r 
Signale sollen auf die S 2 -Mode keinen Einfluss haben, ebenso wenig wie 0.1 
35 Tesla Sinussignale von T 2 auf die Theta-Mode. 

Fiour 23: 



Analog zu Beispiel 2. 
40 Fiaur24: 
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D 1300 (Ubersetzungsoperatoren): 

Falls i und k verschieden sind, sind diese Ubersetzungsoperatoren bspw. fur die 
beiden Slgnale, wie gerade zu Figur 22 dargestellt, gleich Null. 

5 S 2410 (HAC-Prognose): 

Wie vorher. 

Fur die S r Mode: S 2420 (Andern): 

Beispiel: Senden eines Ti 0.5 Tesla 7 Hertz um 3 verschobenen Sinussignal 
1 0 schaltet die SrMode aus. 

Fur die S r Mode: S 2430 (Ermitteln der Sa-Liicke): 

Das Prognose-HAC fur S 2 ergibt sich, da die T r Signale keine direkte Wirkung auf 
die S 2 -Mode haben, aus deren unveranderter Fortsetzung, dieses soli mit dem 
15 S2-lst-HAC ubereinstimmen. 

Fur die SrMode: S 2420 (Andern): 

Beispiel: Frequenzverschieben der S 2 -Mode auf eine mit 7 Hertz 
inkommensurable Grundfrequenz, beispielsweise durch zeitverzogerte 
20 Ruckkopplung des umgepolten Messsignal von S 2 durch T 2 . 

Fiir die S 2 -Mode: S 2430 (Ermitteln der S r Lucke): 

Das Prognose-HAC fiir Si ergibt sich, da die TV-Signale keine direkte Wirkung auf 
die SrMode haben, aus deren unveranderter Fortsetzung, dieses soil mit dem 
25 Srlst-HAC ubereinstimmen. Hier ergibt sich beispielsweise eine Veranderung 
der in Figur 19 (iii) dargestellten Signalkomponente in Richtung Figur 19 (ii). 

S 2440 (Ruckrechnen auf liickenerklarende Funktionalmatrizen): 

Keine Liicke bei SrVeranderung, daher ist die Zelle (1,2) gleich Null. (Mode 1:= 
30 die dargestellte S r Mode, Mode 2:= die dargestellte SrMode). Mehrdeutige 
Liicke bei Zelle (2,1). 

E 400 (Eindeutig?): 

Prufen, ob mehr als eine Funktionalmatrix sich als Ergebnis von S 2440 ergibt. 

35 Wenn man sich in der jeweiligen Anwendung nicht fiir beliebige Zusammenhange 
zwischen Hirnaktivitaten interessiert, sondem nur fQr bestimmte 
Zusammenhange, beispielsweise nur zwischen den Hirnaktivitaten i und j, so 
geniigt es, die (i.j)-te und die G,i)-Zelle der Funktionalmatrix eindeutig zu fullen. Im 
vorliegenden Beispiel interessiert lediglich, ob HAC-erhaltende Veranderungen 

40 der Mode 1 sich ausbreiten, d.h. Mode 2 verandern, hierfur geniigt die 
Betrachtung der Zelle (1,2), und diese hat den Wert Null. 
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S 2221 (Signal andern): 

Beispielsweise Erzeugen von Frequenzverschiebungen, Phasenverschiebungen, 
Amplitudenanderungen, Stationaritat andern, Rauschen hinzufugen, Ljapunov- 
5 Exponenten verandern, u.v.m., was im vorliegenden Beispiel nicht weiter benotigt 
wird. 



D 2410 (Funktionalmatrix): 

Das Ergebnis umfasst jeweils eine Funktionalmatrix, die hier bspw. folgende 



10 Gestalt hat: 



K unbestimmt 



Figur 25: 



Analog zu Beispiel 2. Zusatzlich kann dank der translokalen Kalibrierung 
15 festgestellt werden, dass beispielsweise eine Ausbreitungseinschrankung 
(D4300) der Gestalt: "Durch Ti hervorgerufene Anderungen sollen sich nicht auf 
andere als von Si erfasste Hirnaktivitaten auswirken" erfullt wird. Im vorliegenden 
Beispiel wird das FmB bei beispielsweise auf Null gesunkenem Driving (s. Figur 
19 i) mit Hilfe des o.g. unterbrochenen Sinussignals durch Ti nahe seiner 
20 naturlichen Leistungsspektren ( M volles Driving", s. Figur 19, (ii) und (Hi)) 
stabilisiert, ohne die von S2 erfasste Hirnaktivitat zu storen. Das Ergebnis 
(gemessen von Si) ist in Figur 19 (iv) dargestellt Das Hirnaktivitatsziel der 
Wiederherstellung von Atti, ohne Att 2 zu storen, ist damit erreicht, bei 
Anwendbarkeit des Ishihara-Yoshii-Modell somit auch das Verhaltensziel. 

25 

Erlauterungen und Anmerkungen: 



Zu D 1000: 

Ein Hirnaktivitatsmodell ist im einfachsten Fall eine Menge von Differenzen- 
30 und/oder Differentialgleichungen bspw. der Gestalt 

*(0 = f(x,tJokaleParameter,endogenerInput 9 exogenerInput) 

fur eine Hirnaktivitat, x Hirnaktivitatscharakteristikum, t Zeit. M Exogener Input" 
bezieht sich auf kiinstlich vermoge extrakranial generierter elektromagnetischer 
Felder hervorgerufene Werte von Inputvariablen in den Gleichungen des 

35 Hirnaktivitatsmodells ( im Gegensatz zu ,,endogenem Input", der im Regelfall uber 
normale Sinnes- oder Nervenkanale zustande kommt, beispielsweise durch 
photische Stimulation, durch thalamisches Pacemaking, durch Horen einer 
Symphonie, u.v.m). Exogener Input = U(Transmittersignal), 0 
Ubersetzungsoperator. Die Funktion f richtet sich nach dem verwendeten 

40 Hirnaktivitatsmodell, und bestimmt die Ableitung von x unter Verwendung von x, t, 
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sowie nichtbeobachtbarem endogenem und nichtbeobachtbarem exogenem Input 
fur einen nichtbeobachtbaren Parametersatz. 

Allgemeine Hirnaktivitatsmodelle verwenden statt punktweisem x Pfade und/oder 
5 stochastische Differentialgleichungen unter Einbeziehung von Rauschen. 

Verschiedene Hirnaktivitatsmodelle fiihren zu verschiedenen 
Kalibrierungsergebnissen, mit eventuell verschiedenen resultierenden 
Modifikationen. Fur geeignete Hirnaktivitatsmodelle wird das Hirnaktivitatsziel 
10 zuverlassig erreicht und bei geeignetem Verhaltensmodell das Verhaltensziel. 
Hirnaktivitats- und Verhaltensmodeile sind auch in integrierter Form darstellbar. 

Zu S 2000: 

Die Nichtlinearitat des Systems gestattet es, durch geeignete Variation eines 
15 Transmittersignals moglichst kleiner Feldstarke („Testsignal") die 
nichtbeobachtbaren GroBen inklusive des Obersetzungsoperators 0 zu 
bestimmen. 

Die hier vorgestellte Dynamische Kalibrierung unterscheidet sich sowohl von 
20 vorhandenen mathematischen Methoden zur Parameterschatzung (die in der 
Regel Verteilungsannahmen der als Zufallsvariable betrachteten Parameter 
voraussetzen, und durch Beobachten des Verlaufs von Trajektorien uber die 
Annahme einer zunehmenden Annaherung von aus den Beobachtungen 
errechneten und den wahren Parametern auf die Werte letzterer riickschlieften, 
25 siehe z.B. [7]), als auch von aus den Ingenieurwissenschaften bekannten 
adaptiven Regelungen / Parameterschatzungen (diese sind im wesentlichen nur 
fur lineare zeitinvariante Systeme brauchbar, allgemein tastet man sich als 
Bestandteil der Regelung mit sehr einfachen Funktionen als Stellgrdfien iiber 
eine Minimierung von jdeale Regelgrofce minus tatsachliche Regelgrdfie" voran, 
30 siehe z.B. [8]). 

Dynamische Kalibrierung basiert im Unterschied zur ..passiven Messung" 
(konventionellen nichtinvasiven in-vivo Beobachtung) der dargestellten 
mathematischen Verfahren auf „aktiver Messung" (Obertragung von Signalen in 
das betreffende System, sowie Messung des Systems in Interaction mit dem 

35 Signal). Das Ergebnis aktiver Messung ist vielfach in nichtlinearen 
deterministischen und/oder nichtlinearen stochastischen Systemen (wie dem 
menschlichen Gehirn) deutlich verschieden von den in der Messung linearer 
Systeme iiblichen Impuls- oder Sprungantworten, oder ahnlichem (z.B. 
Stochastische Resonanz). 

40 Im Unterschied zu bekannten adaptiven Regelungen / Parameterschatzungen ist 
die Dynamische Kalibrierung erstens ein im Normalfali von einer Regelung zu 
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separierender und dieser vorangehender Prozess, zweitens liegt bei der 
Dynamischen Kalibrierung das Hauptaugenmerk auf den verwendeten 
Testfunktionen, deren Vielfalt, und damit Trennfahigkeit zwischen verschiedenen 
Satzen von Nichtobservablen der Komplexitat des zu untersuchenden Systems 
5 angemessen ist (im Gegensatz zu einfachen Testfunktionen, wie Heaviside- 
Funktionen o.a., deren Systemantwort vielfach uber lange Zeitraume beobachtet 
werden muss, was im in-vivo Fall ebenfalls nicht zumutbar ist), drittens ist die Art 
der Approximation an einen die Messwerte erklarenden Parametersatz 
grundverschieden (sukzessive Verkleinerung von Teilmengen einer 
10 diskretisierten Parametermannigfaltigkeit, im Gegensatz zum ublichen Soll-lst- 
Vergleich von Reglergrolien). 

Des weiteren liefert Dynamische Kalibrierung zusatzlich zur Feststellung der 
Nichtobservablen die im Falle heterogener und/oder komplexer Systeme 
entscheidenden unbekannten Ubersetzungsoperatoren von Kenngrofien 
15 systemextern erzeugter Felder (z.B. Feldstarke und/oder Frequenzmuster 
extrakranialer Magnetfelder) in exogenen Input der Systemgleichungen. 

ZuD2410: 

Eine Funktionalmatrix ergibt sich als mxm Matrix aller m beobachteten 
20 verschiedenen Hirnaktivitaten, wobei die Matrixzelle my = 0 ist, falls keine 
funktionale Verbindung von der j-ten Aktivitat zur i-ten Aktivrtat zu ermitteln ist. 
Funktionale Verbindungen zwischen Hirnaktivitaten setzen physiologische 
Verbindungen der unterliegenden Himbestandteile voraus, letztere Verbindungen 
gestatten aber keine Aussage uber die Kopplung von Hirnaktivitaten. Mogliche 
25 Funktionalmatrizen sind solche, die nach Parameterdiskretisierung mit den 
Beobachtungen kompatible Losungen fur die Gesamtheit der beobachteten 
Hirnaktivitaten iiefem, wobei in den Gleichungen fur eine Hirnaktivitat zi deren 
endogener Input ersetzt wird durch bspw. Residualinput + Summe aller 
vorgelagerten Aktivitaten zj * jeweilige Einflussstarke* Zj(t-Verzogerung). (,,A 
30 vorgelagert B" heifct „Veranderung von A verandert B", aquivalent hierzu: B 
nachgelagert A.) 

Die bevorzugte Ausfiihrungsform des Verfahren bei lokaler Modifikation umfasst 
folgende Elemente (unter der Annahme, dass ein Cluster von benachbarten 
35 Wilson-Cowan Oszillatoren durch ein Paar von Wilson-Cowan-Gleichungen 
modelliert werden kann): 

a) Datenanalyse mit Fourieranalyse, alternativ auch in Kombination mit 
Waveletanalyse (geeignete Linearkombinationen von Losungen x(t) und y(t) der 
Wilson-Cowan-Gleichungen als Basis) (Grundlagen siehe [13]), 
40 Phasenraumeinbettung, sowie statistischen Verfahren, urn bezuglich der o.g. 
Analysemethoden die Wahrscheinlichkeit der Ubereinstimmung von 
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Komponenten des gemessenen Signals mit einer Aquivalenzklasse von 
Losungen der Wilson-Cowan-Gleichungen festzustellen. 
Testsignale vorzugsweise linear kombiniert aus zeitverzogerter Riickkopplung 
von Messsignalen und/oder Losungen der Wilson-Cowan-Gleichungen, d.h. 

5 beispw. M ' ~ *) - y(t - t)) . 

b) Lokale Kalibrierung, 

c) Modifikation wie beschrieben, 

d) laufende Anfallswamung nebst automatischer Gegensteuerung/-regelung wie 
beispw. in [9]. 

10 

Analog ergibt sich die bevorzugte Ausfiihrungsform fur translokale Phanomene, 
wie oben dargestellt, plus translokale Kalibrierung. 

Falls keine Kalibrierung bendtigt wird und die Limit-Cycle-Oszillatoren nur 
15 schwach gekoppelt sind und die Annahme gilt, dass sich Cluster von schwach 
gekoppelten Limit-Cycle-Oszillatoren durch Cluster von Phasenoszillatoren auch 
fur kurze Zeitraume ersetzen lassen, konnen statt Limit-Cycle-Oszillatoren 
(Losungen der Wilson-Cowan-Gleichungen) gekoppelte Phasenoszillatoren 
(siehe beispielsweise [23]) verwendet werden. Interventionen basieren auf 
20 Phasenreset und anschlieliender extemer Versklavung nebst laufender Messung 
eines Signals, sodass der Cluster sich zielgemalJ entwickelt. Diese 
Ausfiihrungsform ist vorzugsweise bei eng umschriebenen Hirnarealen 
anzuwenden, damit der Reset keine unerwunschten Storungen hervorruft. 

25 Fur einen Messkanal wird die Phasenraumeinbettung illustriert: 

Fig. 26 zeigt die EEG-Daten eines Kanals, wobei die x-Achse Zeiteinheiten 
(1/128 eines Sekunde) darstellt, die y-Achse Spannungsdifferenzen (zwischen 
der Mess- und der Referenzelektrode zum jeweiligen Zeitpunkt). Im einfachsten 
Fall erfolgt Feststellung des Oder der Orbits beispielsweise gemafc [29]. Allgemein 

30 erfolgt die Einbettung uber Delay-Koordinaten (siehe [27]), die Delay-Time ist, wie 
dort dargestellt, beispielsweise als erster Nulldurchgang der 
Autokorrelationsfunktion zu ermitteln. Einbettung eines Datenfensters fur t = 130 
und T = 86 (s. Figur 27). Hierbei gibt die x-Achse fur ein mit t endendes 
Fenster der Lange 32 ax(t-86)+p-y(t-86) wieder, und die y-Achse a x(t)+B y(t). In 

35 Fig. 27 ist alpha auf 1 und beta auf -1 gesetzt. Da Mehrfachbelegungen 
desselben Punktes nicht dargestellt werden, ist die Anzahl der Punkte in Figur 27 
kleiner als 32. Mit gangigen Verfahren (Bayes, Minimax, u.v.a) wird nun nach 
Reskalierung die Obereinstimmung der empirischen Phasenraumbelegung mit 
jeder (auf dem endlichen Parametergitter) der Wilson-Cowan-Losungen fur jeden 

40 (auf den diskretisierten Einheitintervallen) Mischungswinkel ermittelt. Figur 28 
stellt mit gleichen Achsenbezeichnungen wie Figur 27 und gleichen 
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Mischungswinkeln wie oben die Losungen fur tau=10, a=22, b=20, c=14, d=3, 
rho x =1.5 bei verschiedenem rho y dar Rauten fur rho y =-1.5, Quadrate fiir y=0. 
Figur 29 zeigt die Bedeutung der Mischungswinkel: hier sind bei gleichen 
Achsenbezeichnungen wie vorher die Losungen fiir tau=10, a=22, b=20, c=14, 
5 d=3, rhOx=1.5 , rho y =0 mit alpha=1 und beta=-1 (Raute) sowie beta=-0.7 
(Quadrate), und beta=-0.3 (Kreuze) gegenubergestellt. 

In einer vorteilhaften Ausgestaltung des Verfahrens werden statt der Wilson- 
Cowan-Gleichungen stochastische Wilson-Cowan-Gleichungen verwendet, das 
1 0 heiftt, Px = Ac.det + Rauschen , p y = p yM + Rauschen . nde t" steht fur 

den deterministischen Inputanteil, Rauschen fiir den stochastischen Input, das 
gleiche Verfahren wie vomer wird angewandt. Zusatzlich erfolgt die Feststeilung 
der Rauschlevels. Sinnvoll ist hier ebenfalls die Verwendung von verrauschten 
Testsignalen und Messung von Signal-To-Noise Ratios. 

15 

In einer weiteren vorteilhaften Ausgestaltung des Verfahrens werden in den 
Wilson-Cowan-Gleichungen fur x und y unterschiedliche Zeitkonstanten 
? verwendet und / oder unterschiedliche Multiplikatoren M . Deren Kalibrierung 
erfolgt wie bereits fiir die anderen Parameter beschrieben. 

20 

In einer weiteren vorteilhaften Ausgestaltung des Verfahrens werden in den 
Wilson-Cowan-Gleichungen jedes Wilson-Cowan-Neurons eines Clusters die 
anderen Neuronen des Clusters iiber Mean-Field-Coupling eingekoppelt, das 
heiRtfurallei = 1..N 

N N 

25 tx, = -x, +S(ax i -by, +p x>i +s £ Xj + 5 £ yj ) sowje 

N N 

Ty i =-y i +S(cx i -dy i +p yJ +T J £ xj+v X yj) mit 

neuen, zu ermittelnden Nichtobservablen. Analog ist dies fur Cluster 
von stochastischen Wilson-Cowan-Neuronen, analog mit oder ohne 
zeitverzogerten Input ausfiihrbar. 

30 

In einer weiteren vorteilhaften Ausgestaltung des Verfahrens wird die Gestalt der 
sigmoiden Funktion der Neuronen eines Clusters mit Hilfe von verrauschten 
Testfunktionen kalibriert, da der Cluster je nach sigmoider Funktion Signal bzw. 
Rauschen verschieden verstarkt bzw. unterdrOckt (siehe fiir einfache Netzwerke 
35 beispw. [28]). 
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In einer weiteren vorteilhaften Ausgestaltung des Verfahrens werden 
Testfunktionen nacheinander jeweils umgekehrter Polung verwendet, das heifit 
zunachst Verwendung von eines Transmittersignals xi(t), danach von -xi(t). 

5 In einer vorteilhaften Ausgestaltung wird die Kalibrierung nach n > 0 Testsignalen 
abgebrochen („Schnellkalibrierung"), ein mit den Messwerten kompatibler Satz 
von Parametern sowie endogenen und exogenen Inputs als 
Kalibrierungsergebnis vorlaufig festgelegt. Sich im Rahmen der Modifikation 
(S3000) ergebende zusatzliche Messwerte werden laufend zur Anpassung des 
1 0 oben genannten Satzes verwendet. 

In einer vorteilhaften Ausgestaltung werden Einflussannahmen und -invarianten 
(D1049, D1050)validiert. 

15 In einer vorteilhaften Ausgestaltung wird das EM Verfahren mit sensorischem 
Input (akustisch, optisch, u.a.) kombiniert. Beispiel: Instruktion an den Anwender 
in Beispiel 3, mit Beginn der Kalibrierung eine mentale Aufgabe (s. beispw. [20]) 
durchzufuhren, wodurch das Auffinden geeigneter Fmd -Moden erleichtert wird. 
Anderes Beispiel: Akustisches Pacemaking zur translokalen Kalibrierung 

20 auditorischer Hirnprozesse. 

In einer vorteilhaften Ausgestaltung wird der Effekt des Verfahrens mit Hiife 
geeigneter, fur den einzelnen Anwender zu per externer Validierung zu 
ermittelnder Wiederholungsraten uber die Dauer der einzelnen Anwendung 
25 hinaus stabiiisiert. Diese Auswirkung geeigneter Wiederholungsraten ist bspw. 
fur Anwendungen von TMS (transkranialer Magnetstimulation mit Magnetfeldern 
von 1-2 Tesla) bekannt (s. beispw. [3]). 

In einer vorteilhaften Ausgestaltung wird zusammen mit der Kalibrierung 
30 Sensoroptimierung betrieben: Zuschalten vorher inaktiver Sensoren und/oder 
Positionsveranderung von Sensoren und/oder Ausrichtungsveranderung von 
Sensoren solcherart, dass diese identifizierte und/oder nichtidentifizierte 
Hirnaktivitaten vom betreffenden, ggf. neu ausgerichteten Sensor besonders gut 
erfasst werden. Zur Sensoroptimierung sind herkommliche 
35 Optimierungsverfahren geeignet. 

In einer vorteilhaften Ausgestaltung wird zusammen mit der Kalibrierung 
Sensoroptimierung durch Zusammenschalten existierender Sensoren zu 
virtuellen Sensoren betrieben, solcherart dass Hirnaktivitaten vom betreffenden 
40 virtuellen Sensor besonders gut erfasst werden. Beispiel: Si, S2, S3, S4, uber 
eine Funktion f(Si, S 2 , S 3 , S 4 )=0.3*Si+0.5*S2+0.02*(S3) 2 +0.5*sin(Si+S 4 ) zu Svirtueii 
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zusammengeschaltet. Ein geeignetes f wird fur jede Teilmenge von Sensoren mlt 
Hilfe herkommliche Optimierungsverfahren ermittelt 

In einer vorteilhaften Ausgestaltung wird zusammen mit der Kalibrierung 
5 Transmitteroptimierung betrieben: Zuschalten vorher inaktiver Transmitter 
und/oder Positionsveranderung von Transmittern und/oder 
Ausrichtungsveranderung von Transmittern solcherart, dass diese Himaktivitaten 
vom betreffenden, ggf. neu ausgerichteten Transmitter besonders gut beeinflusst 
werden. Zur Transmitteroptimierung sind herkommliche Optimierungsverfahren 
10 geeignet. 

In einer vorteilhaften Ausgestaltung wird zusammen mit der Kalibrierung 
Transmitteroptimierung durch Zusammenschalten existierender Transmitter zu 
virtuellen Transmittern betrieben, solcherart dass Himaktivitaten vom 
15 betreffenden virtuellen Transmitter besonders gut beeinflusst werden. Das 
Optimum (der virtuelle Transmitter) ergibt sich aus den Ubersetzungsoperatoren 
mit Hilfe herkommiicher Optimierungsverfahren 

In einer vorteilhaften Ausgestaltung von S 2130 wird die Dekomposition fur n 
20 Sensoren (n ganze Zahl zwischen 1 und der Maximalzahl der verwendeten 
Sensoren) gleichzeitig durchgefuhrt, d.h. die zu zerlegende Zeitreihe ist ein n- 
Tupel von Ein-Sensor-Zeitreihen. 

In einer vorteilhaften Ausgestaltung von S 2130 wird die Dekomposition fur ein n- 
25 Tupel von Ein-Sensor-Zeitreihen mit Pursuit-Verfahren durchgefuhrt (z.B. 
Matching Pursuit, s.[13]). 

In einer weiteren vorteilhaften Ausgestaltung von S 2130 erfolgt die 
spatiotemporale Dekomposition des n-Tupels von Ein-Sensor-Zeitreihen 
30 beispielsweise mit dem Karhunen-Loeve-Verfahren (s. z.B. [17], [18]) und/oder 
Independent Component Analysis („ICA", s. bspw. [22]). 

In einer weiteren vorteilhaften Ausgestaltung von S 2130 erfolgt die 
Dekomposition nach Einbettung in einen Metaphasenraum, wodurch stationare 
35 und nichtstationare Anteile getrennt werden (s. beispw. [1 5]). 

In einer vorteilhaften Ausgestaltung von S 2130 erfolgt die Dekomposition parallel 
nach mehreren Verfahren, wobei nur die (mit einer a priori definierten 
Gewichtung) von einer Mehrzahl dieser Verfahren identifizierten Moden zur 
40 Weiterverarbeitung (u.a. Relevanzpriifung) akzeptiert werden. 
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In einer vorteilhaften Ausgestaltung von D 2000 wird die Veranderung oder 
Nichtveranderung von Wahmehmung und/oder Wahmehmungsfahigkeit und/oder 
Wahrnehmungsbereitschaft als Verhaltensziel gesetzt. Beispiel hierfur ist u.a. die 
Veranderung der Diskriminationsfahigkeit fur bestimmte Reize oder Klassen von 
5 Reizen. 

In einer vorteilhaften Ausgestaltung von D 2000 wird die Veranderung oder 
Nichtveranderung von Handlungen und/oder Handlungsfahigkeit und/oder 
Handlungsbereitschaft als Verhaltensziel gesetzt. Beispiel hierfur ist u.a. die 
10 Verbesserung von Reaktionsgeschwindigkeit. 

In einer vorteilhaften Ausgestaltung von D 2000 wird die Veranderung oder 
Nichtveranderung von Aktivierung und/oder Aktivierungsfahigkeit als 
Verhaltensziel gesetzt. 

15 

In einer vorteilhaften Ausgestaltung von D 2000 wird die Veranderung oder 
Nichtveranderung von Motivation und/oder Motivierungsfahigkeit und/oder 
Motivierungsbereitschaft als Verhaltensziel gesetzt. 

20 In einer vorteilhaften Ausgestaltung von D 2000 wird die Veranderung oder 
Nichtveranderung von Aufmerksamkeit und/oder Aufmerksamkeitsfahigkeit als 
Verhaltensziel gesetzt. 

In einer vorteilhaften Ausgestaltung von D 2000 wird die Veranderung oder 
25 Nichtveranderung von Gedachtnis und/oder Gedachtnisinhalten und/oder 
Gedachtniszugriff als Verhaltensziel gesetzt. 

In einer vorteilhaften Ausgestaltung von D 2000 wird die Veranderung oder 
Nichtveranderung von Lemen und/oder Lemfahigkeit und/oder Lernbereitschaft 
30 als Verhaltensziel gesetzt. 

In einer vorteilhaften Ausgestaltung von D 2000 wird die Veranderung oder 
Nichtveranderung von Bewusstsein als Verhaltensziel gesetzt. 

35 In einer vorteilhaften Ausgestaltung von D 2000 wird die Veranderung oder 
Nichtveranderung von Emotionen und/oder Emotionsfahigkeit und/oder 
Emotionsbereitschaft als Verhaltensziel gesetzt. 

In einer vorteilhaften Ausgestaltung von D 2000 wird die Veranderung oder 
40 Nichtveranderung von Appetenzen und/oder Aversionen als Verhaltensziel 
gesetzt. 
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In einer vorteilhaften Ausgestaltung von D 2000 wird die Veranderung Oder 
Nichtveranderung von Denken und/oder Denkfahigkeit und/oder Denkbereitschaft 
als Verhaltensziel gesetzt. („Denken" wird hier als „Durchfuhrung mentaler 
5 Prozesse" definiert.). 

In einer vorteilhaften Ausgestaltung von D 2000 wird die Veranderung oder 
Nichtveranderung von Verhaltensablaufen als Verhaltensziel gesetzt. 

10 In einer vorteilhaften Ausgestaltung von D 2000 wird die Veranderung oder 
Nichtveranderung von Verhaltenskorrelationen als Verhaltensziel gesetzt. 

In einer vorteilhaften Ausgestaltung von D 2000 werden mehrere kompatible 
Verhaltensziele kombiniert. 

15 

In einer vorteilhaften Ausgestaltung von D 2000 werden mehrere nicht 
notwendigerweise kompatible Verhaltensziele hierarchisch kombiniert, d.h. es 
wird eine Prioritatsliste der Verhaltensziele aufgestellt und bei gegebenenfalls 
einander widersprechenden Verfahrenschritten die zum niedriger priorisierten 
20 Verhaltensziel gehorigen Verfahrensschritte spater oder gar nicht ausgefuhrt. 

In einer vorteilhaften Ausgestaltung von D 2000 wird zu anderen Verhaltenszielen 
dasjenige der Anfalisfreiheit mit oberster Prioritat hinzugefiigt, was eine 
Anwendung von [9] auf nicht medizinisch vorbelastete Anwender beinhaltet. 

25 

In einer vorteilhaften Ausgestaltung von D 3000 wird das benotigte 
Verhaltensmodell der aulierst umfangreichen Literatur entnommen (z.B. fur die 
vorgestellten Beispiele aus [11], [20], [21] jeweils Dutzende von 
Veroffentlichungen mit gleichem/ahnlichem Gegenstand). 

30 

In einer vorteilhaften Ausgestaltung von D 3000 werden mehrere 
Verhaltensmodelle parallel verwendet, wobei nur diejenigen Steuerungs- oder 
Regelungsschritte (S 3000) durchgefuhrt werden, die (mit einer a priori definierten 
Gewichtung) mit einer Mehrzahl dieser Modelle kompatibel sind 

35 

In einer vorteilhaften Ausgestaltung erfolgt mit Hilfe von i.d.R. psychologischen 
Tests externe Validierung der Verhaltensanderung/-aufrechterhaltung fur den 
einzelnen Anwender wahrend und/oder nach der EM Anwendung. Hierdurch 
wird gepKift, in welchem Umfang das verwendete, i.d.R. statistisch validierte 
40 Verhaltensmodell auf den betreffenden Anwender zutrifft. 
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In einer vorteilhaften Ausgestaltung von D 1000 werden mehrere generische 
Hirnaktivitatsmodelle parallel verwendet, wobei nur diejenigen Steuerungs- oder 
Regelungsschritte (S 3000) durchgefuhrt werden, die (mit einer a priori definierten 
Gewichtung) mit einer Mehrzahl dieser Modelle kompatibel sind. 

5 

In einer vorteilhaften Ausgestaltung des Verfahrens werden in den jeweiligen 
Schritten einzelne Sensoren und/oder einzelne Transmitter durch Mengen von 
Sensoren bzw. Mengen von Transmittem ersetzt. Dieses dient zur Feststellung 
und kontrollierten Modifikation nichtlokaler spatiotemporaler Moden, wie sie bspw. 
1 0 bei Jirsa-Haken-Modellen auftreten (s. bspw. [1 8]). 

In einer vorteilhaften Ausgestaltung von D 4000 werden die minimalen 
Apparaturvoraussetzungen iiberschritten, urn Zusatzinformationen erhalten zu 
konnen, die bei spateren Anwendungen mit moglicherweise geanderten 
15 Verhaltenszielen von Nutzen sein konnen (z.B.: grundsatzliche Anwendung des 
10/20 Oberflachen EEGs, s. z.B. [16]). Der Vorteil dieser Ausgestaltung liegt in 
moglicher Standardisierung von Tellen des Apparates sowie teilweise moglicher 
Verwendung handelsublicher Hardware. 

20 In einer weiteren vorteilhaften Ausgestaltung von D 4000 werden die minimalen 
Apparaturvoraussetzungen uberschritten, urn zusatzliche Steuerungs- oder 
Regelungsmdglichkeiten zu haben ( z.B.: zu jeder Nicht-Referenzelektrode des 
10/20 Oberflachen EEGs eine unmitteibar frontal vor dieser Elektrode befindliche 
Spule). Der Vorteil dieser Ausgestaltung liegt in moglicher Standardisierung von 

25 Teilen des Apparates. 

In einer vorteilhaften Ausgestaltung von S 3000 werden als die Kontrollvariable 
eingesetzten exogenen EM Felder durch sensorische Inputs und/oder 
biochemisch wirksame Substanzen erganzt, wobei sich an den dargestellten 
30 Verfahrensschritten grundsatzlich nichts andert. 

In weiteren vorteilhaften Ausgestaltungen des Verfahrens werden als Test- und 
Transmittersignale Signale mit stochastischen Anteilen verwendet (insbesondere 
bei verwendeten Hirnaktivitatsmodellen mit expliziten stochastischen 
35 Komponenten, wie bspw. [25], einer Erweiterung von [18], mit dem Ziel der 
Nutzbarmachung von stochastischer Resonanz und Coherence Resonance, 
siehe bspw. [26] fur das FitzHugh-Nagumo-Modell). 

In weiteren vorteilhaften Ausgestaltungen des Verfahrens werden als Test- und 
40 Transmittersignale Signale mit chaotischen Anteilen verwendet (insbesondere bei 
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verwendeten Hirnaktivitatsmodellen mit expliziten chaotischen Komponenten, wie 
bspw. Erweiterungen von [1] mit lokal stark gekoppelten Oszillatoren). 

In einer weiteren vorteilhaften Ausgestaltung des Verfahrens werden inaktive 
5 Moden durch Transmittersignale direkt aktiviert, d.h. durch den unmittelbaren 
Effekt der Signale auf die Mode. 

In einer weiteren vorteilhaften Ausgestaltung des Verfahrens werden Moden 
durch Transmittersignale indirekt beeinflusst, d.h. uber direkten Einfluss auf 
10 Moden, was wiederum (s. Funktlonalmatrix) Auswirkungen auf die indirekt zu 
beeinflussenden Moden hat. 

In einer weiteren vorteilhaften Ausgestaltung des Verfahrens werden 
Bifurkationspunkte kontrolliert, d.h. es werden qualitative Veranderungen EM 
15 Himaktivitat vorgenommen (wie etwa in Beispiel 3, wo der ohne Driving inaktive 
neurale Oszillator eingeschaltet wurde). 

In einer weiteren vorteilhaften Ausgestaltung des Verfahrens werden inaktive 
Moden aktiviert. 

20 

In einer weiteren vorteilhaften Ausgestaltung des Verfahrens erfolgt 
Amplitudenverstarkung durch Aktivierung inaktiver Moden und/oder 
Synchronisierung von Moden, analog Abschwachung durch Ausschaltung aktiver 
Moden und/oder Desynchronisierung von Moden. 

25 

In einer weiteren vorteilhaften Ausgestaltung des Verfahrens erfolgt 
alterskorrelierter Verflachung (fur EEG s. beispielsweise [12]) entgegenwirkende 
Amplitudenverstarkung. 

30 In einer weiteren vorteilhaften Ausgestaltung des Verfahrens erfolgt die 
Dekomposition nichtstationarer Signale zunachst in einen stationaren und einen 
nichtstationaren Teil (s. bspw. [15]), wobei die verwendeten Test- und 
Transmittersignale dann mit dem nichtstationaren Anteil zu modulieren sind. 

35 In einer weiteren vorteilhaften Ausgestaltung des Verfahrens erfolgt die 
Identifikation von Anwendern auf der Grundlage zeitstabiler Kalibrierungsdaten. 
Bereits herkdmmllch spektralzerlegter EEG-Daten haben fur ein Individuum 
zeitstabile Komponenten (nachgewiesen bis zu funf Jahren Dauer bspw. [10]). 
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Anspruche: 

1. Verfahren zur elektromagnetischen Modifikation von Hirnaktivitat zum 
Erreichen eines vorgegebenen Verhaltenszieles, unter Verwendung eines den 

5 Zusammenhang zwischen Verhalten und den Verlauf von 

Hirnaktivitatscharakteristika beschreibenden Verhaltensmodells, sowie eines 
Hirnaktivitatsmodells, welches Hirnaktivitat quantitativ beschreibt, wobei 
Hirnaktivitatscharakteristika und deren Verlaufe aus der quantitativen 
Beschreibung der Hirnaktivitat ableitbar sind, umfassend 
1 0 - Ermitteln mittels des Verhaltensmodells, welcher Soil- Verlauf 

bestimmter Hirnaktivitatscharakteristika dem Erreichen des Verhaltenszieles 
entspricht, 

Erzeugung exogener elektrischer und/oder magnetischer Felder, und 
Messung elektromagnetischer Hirnaktivitat nebst Berechnung von 
15 Hirnaktivitatscharakteristika zum Erreichen des Soll-Verlaufs in einem Steuer- 
oder Regelkreis, der Berechnungen mit dem Hirnmodell beinhaltet. 

2, Verfahren nach Anspruch 1 , 

wobei ein generisches Hirnaktivitatsmodell verwendet wird, das durch Ermitteln 
20 von Nichtobservablen des generischen Hirnaktivitatsmodells fur einen 

bestimmten Anwenderzu einem spezifischen Hirnaktivitatsmodell kalibriert wird, 
und der Einfluss exogener elektrischer und/oder magnetischer Felder auf 
Hirnaktivitat fur diesen Anwender ermittelt wird. 

25 3. Verfahren nach Anspruch 1 oder 2, 

wobei als Hirnaktivitatscharakteristika die Hirnaktivitat beschreibende 
elektromagnetische Grolien oder aus elektromagnetischen Groften abgeleitete 
Hilfsgrofien verwendet werden. 

30 4. Verfahren nach Anspruch 3, 

wobei die elektromagnetischen GrofJen extra- und/oder intrakranial gemessene 
Potenziale und/oder Strome und/oder Magnetfelder umfassen. 

5. Verfahren nach Anspruch 3 oder 4, 

35 wobei die Hirnaktivitatscharakteristika Frequenzen elektromagnetischer Grolien 
und/oder Koeffizienten bzgl. einer Fourier- oder Wavelet-Basisentwicklung oder 
Karhunen-Loeve-Koeffizienten von Zeitreihen der elektromagnetischen Groften 
und/oder stochastische Kenngroften, wie z.B. Stationaritat, Momente, 
Korrelationen, und/oder nichtlineare KenngroUen, wie z.B. Stabilitat Ljapunov- 

40 Exponenten, Grenzzyklen, minimale Phsenraumdimension, umfassen. 
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6. Verfahren nach einem der Anspruche 1 bis 5, 

wobei zum Berechnen des Soll-Verlaufs aus dem Verhaltensziel ein 
Verhaltensmodell verwendet wird, wie z.B. Varianten des Davidson-Modells, das 
Ishihara-Yoshie-Modell oderdas Chaos-Anfallsmodell gemaB [9], 

5 

7. Verfahren nach einem der Anspruche 2 bis 6, 

wobei die Nichtobservablen fur einen bestimmten Anwender kalibriert werden. 

8. Verfahren nach einem der Anspruche 1 bis 7, 

10 wobei tier Einfluss eines exogenen elektrischen und/oder magnetischen Feldes 
auf die Hirnaktivitat eines bestimmten Anwenders ermittelt wird. 

9. Verfahren nach einem der Anspruche 2 bis 8, 

wobei das Berechnen des Einflusses des exogenen elektrischen und/oder 
15 magnetischen Feldes und das Ermitteln der Nichtobservablen gekoppelt 
ausgefuhrt werden. 

10. Verfahren nach einem der Anspruche 2 bis 9, 

wobei beim Kalibrieren zunachst eine Dekomposition (S 2130) von Ist-Verlaufen 
20 von Himaktivitatscharakteristika in HAOModen derart vorgenommen wird, dass 
diesen HAOModen aus dem Himaktivitatsmodell errechnete Verlaufe von 
Himaktivitatscharakteristika entsprechen. 

1 1 . Verfahren nach Anspruch 1 0, 

25 wobei anhand dieser Moden Satze von Parametern und endogenen Inputs 

und/oder exogenen Inputs ermittelt werden, die gemafi dem Himaktivitatsmodell 
der jeweiligen Mode zugrunde liegen konnen. 

12. Verfahren nach einem der Anspruche 1 bis 1 1 , 

30 wobei ein Testsignal berechnet wird (S 2160), wobei das Testsignal ein exogenes 
elektrisches und/oder magnetisches Feld ist, das an den Anwender angelegt 
wird, wobei dessen Hirnaktivitat gemessen wird und die nach oder bei Anlegen 
des Testsignals gemessene Hirnaktivitat zusammen mit vorheriger Hirnaktivitat 
ausgewertet wird und anhand des Auswertungsergebnisses mit den Messwerten 

35 kompatible Parametersatze nebst endogenen Inputs nebst Einflussen des 
Testsignals auf die betrachtete Mode ermittelt werden. 

13. Verfahren nach Anspruch 12, 

wobei das Testsignal solange verandert und erneut angelegt wird, bis sowohl 
40 Parametersatz, endogener Input als auch Einfluss des Testsignals auf die 
betrachtete Mode eindeutig feststehen, wodurch die Mode kalibriert ist. 
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14. Verfahren nach einem der Anspruche 12 und 13, 

wobei die Kalibrierung von Himaktivitat durch Kalibrierung von Moden sowie 
Wechselwirkungen zwischen Moden erfolgt, ersteres und oder beides unter 
5 Berucksichtigung exogenen Inputs. 

15. Verfahren nach einem der Anspruche 1 bis 14, 

wobei eine Schnellkalibrierung bereits nach Senden eines oder weniger 
Testsignale einen Satz von mit den Messwerten kompatiblen Parametern, 
10 endogenen Inputs, exogenen Inputs auswahlt. 

16. Verfahren nach Anspruch 15, 

wobei anhand des kalibrierten Hirnaktivitatsmodells berechnet wird, welcher 
Verlauf einer Kontrollvariable geeignet ist, den Ist-Verlauf der 
1 5 Hirnaktivitatscharakteristika in den Soll-Verlauf zu uberfuhren, wobei die 

Kontrollvariable den dem exogenen magnetischen und/oder elektrischen Feld 
entsprechenden exogenen Input im Hirnaktivitatsmodell reprasentiert. 

17. Verfahren nach Anspruch 16, 

20 wobei der zu berechnende Verlauf der Kontrollvariablen durch Anwendung von 
Signalen bekannter Wirkung vereinfacht wird. 

18. Verfahren nach einem der Anspruche 2 bis 17, 

wobei fur mehrere Messbereiche lokale Moden kalibriert werden, 
25 deren Wechselwirkungen durch sukzessives Verandern einzelner lokaler Moden 
ermittelt werden (S 2420, S 2430, S 2440, S 2221 ). 

19. Verfahren nach einem der Anspruche 1 bis 18, 

wobei erganzend zu den elektrischen und/oder magnetischen Feldern 
30 sensorische Inputs und/oder biochemisch wirksame Substanzen eingesetzt 
werden. 

20. Verfahren nach einem der Anspruche 1 bis 19, 

wobei Himaktivitat mit Hilfe einer geeigneten Wiederholungsrate uber die Dauer 
35 der Anwendung hinaus verandert und oder stabilisiert wird. 

21. Verfahren nach einem der Anspruche 1 bis 20, 

wobei zur Messung der elektromagnetischen Himaktivitat mehrere Sensoren 
zusammengeschaltet werden. 

40 
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22. Verfahren nach einem der Anspruche 1 bis 21, 

wobei zum Erzeugen der Felder ein Transmitter verwendet wird, der bzgl. seiner 
Position und Ausrichtung veranderbar ist. 

5 23. Verfahren nach einem der Anspruche 1 bis 22, 

wobei zum Erzeugen der Felder Transmitter verwendet werden, wobei mehrere 
Transmitter zusammen geschaltet werden. 

24. Verfahren nach einem der Anspruche 10 bis 23, 

1 0 wobei Dekomposition einer n-elementigen Zeitreihe in spatiotemporale Moden 
vorgenommen wird. 

25. Verfahren nach einem der Anspruche 10 bis 24, 

wobei die Dekomposition mit Pursuit-Verfahren durchgefuhrt wird. 

15 

26. Verfahren nach einem der Anspruche 10 bis 25, 

wobei die Dekomposition nach transformierten Outputfunktionen des Hirnmodells 
erfolgt, wobei die jeweilige Transformation eine Signalveranderung zwischen den 
an der Hirnaktivitat beteiligten Neuronen und dem betreffenden Sensor darstellt. 

20 

27. Verfahren nach einem der Anspruche 10 bis 26, wobei die Dekomposition 
nach Einbettung in einen Metaphasenraum erfolgt, wodurch stationare und nicht- 
stationare Anteile trennbar werden. 

25 28. Verfahren nach einem der Anspruche 1 0 bis 27, 

wobei die Dekomposition mehrfach mit unterschiedlichen Verfahren ausgefuhrt 
wird, wobei die mehreren Satze von Moden in Abhangigkeit von den Verfahren 
gewichtet werden und nur die Moden fur das weitere Verfahren verwendet 
werden, deren gesamte Gewichtung uber einem bestimmten Schwellenwert liegt. 

30 

29. Verfahren nach einem der Anspruche 1 bis 28, 

wobei das Verhaltensziel eines oder mehrere der nachfolgend angegebenen 

Verhaltensziele umfasst: 

Veranderung oder Nichtveranderung von Wahmehmung und/oder 
35 Wahrnehmungsfahigkeit und/oder Wahrnehmungsbereitschaft, Veranderung oder 
Nichtveranderung der Diskriminationsfahigkeit fur bestimmte Reize oder Klassen 
von Reizen, 

Veranderung oder Nichtveranderung von Handlungen und/oder 
Handlungsfahigkeit und/oder Handlungsbereitschaft, 
40 Veranderung oder Nichtveranderung von Reaktionsgeschwindigkeit, 
Veranderung oder Nichtveranderung von Aktivierung und/oder 
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Aktivierungsfahigkeit, 

Veranderung oder Nichtveranderung von Motivation und/oder 
Motivierungsfahigkeit und/oder Motivierungsbereitschaft, 
Veranderung oder Nichtveranderung von Aufrnerksamkeit und/oder 
5 Aufmerksamkeitsfahigkeit, 

Veranderung oder Nichtveranderung von Gedachtnis und/oder 
Gedachtnisinhalten und/oder Gedachtniszugriff, 

Veranderung oder Nichtveranderung von Lernen und/oder Lernfahigkeit und/oder 
Lembereitschaft, 
1 0 Veranderung oder Nichtveranderung von Bewusstsein, 

Veranderung oder Nichtveranderung von Emotionen und/oder Emotionsfahigkeit 
und/oder Emotionsbereitschaft, 

Veranderung oder Nichtveranderung von Appetenzen und/oder Aversionen, 
Veranderung oder Nichtveranderung von Denken und/oder Denkfahigkeit 
1 5 und/oder Denkbereitschaft, 

Veranderung oder Nichtveranderung von Verhaltensablaufen, 
Veranderung oder Nichtveranderung von Verhaltenskorrelationen. 

30. Verfahren nach Anspruch 30, 

20 wobei mehrere nicht notwendigerweise kompatible Verhaltensziele hierarchisch 
kombiniert werden. 

31 . Verfahren nach Anspruch 30, 

wobei dem Verhaltensziel der Anfallsfreiheit in der Hierarchie eine hochste 
25 Prioritatzugeordnetwird. 

32. Verfahren nach einem der Anspriiche 1 bis 31, 

wobei mehrere gewichtete Verhaltensmodelle parallel verwendet werden, wobei 
nur diejenigen Steuerungs- oder Regelungsschritte (S 3000) durchgefuhrt 
30 werden, die mit einer Mehrzahl dieser Modelle kompatibel sind. 

33. Verfahren nach einem der Anspriiche 1 bis 32, 

wobei eine externe Validierung der VerhaltensanderungAaufrechterhaltung fur 
den einzelnen Anwender durchgefuhrt wird. 

35 

34. Verfahren nach einem der Anspruche 1 bis 32, 

wobei mehrere gewichtete Hirnaktivitatsmodelle parallel verwendet werden, 
wobei nur diejenigen Steuerungs- oder Regelungsschritte (S 3000) durchgefuhrt 
werden, die mit einer Mehrzahl dieser Modelle kompatibel sind. 

40 
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35. Verfahren nach einem der Anspriiche 2 bis 34, 

wobei Kalibrierungsanforderungen fur Einzelsensoren durch Kalibrierungs- 

anforderungen fur Mengen von Sensoren ersetzt werden. 

5 36. Verfahren nach einem der Anspriiche 1 bis 35, 

wobei das Erzeugen der Felder und das Messen alternierend durchgefuhrt 
werden. 

37. Verfahren nach einem der Anspriiche 1 bis 35, 

10 wobei das Erzeugen der Felder und das Messen gleichzeitig durchgefuhrt 
werden. 

38. Verfahren nach einem der Anspriiche 14 bis 37, 

wobei das zu verwendende Testsignal lineare Oder nichtlineare und/oder 
1 5 chaotische und/oder stochastische und/oder nichtstationare Anteile aufweist. 

39. Verfahren nach einem der Anspriiche 17 bis 38, 

wobei die Identifikation von Anwendern auf der Grundlage zeitstabiler 
Kalibrierungsdaten erfolgt. 

20 

40. Verfahren nach einem der Anspriiche 17 bis 39, 

wobei Kalibrierungs- und Modifikationsschritte fQr Gruppen von Sensoren und 
Transmittern sequentiell oder parallel durchgefuhrt werden. 

25 41 . Verfahren nach einem der Anspriiche 1 7 bis 40, 

wobei Bifurkationspunkte kontrolliert werden, und somit unter anderem inaktive 
Moden aktiviert werden. 

42. Verfahren nach einem der Anspriiche 18 bis 41, 

30 wobei Amplitudenverstarkung alterskorrelierter Aktivitatsverflachung erfolgt. 

43. Verfahren nach einem der Anspriiche 1 bis 42, 
wobei es automatisch ausgefiihrt wird. 

35 44. Vorrichtung zum Erreichen eines vorgegebenen Verhaltenszieles, unter 

Verwendung eines den Zusammenhang zwischen Verhalten und den Verlauf von 
Himaktivitatscharakteristika beschreibenden Verhaltensmodells, sowie eines 
Hirnaktivitatsmodells, welches Himaktivitat quantitativ beschreibt, wobei 
Himaktivitatscharakteristika und deren Verlaufe aus der quantitativen 

40 Beschreibung der Himaktivitat ableitbar sind, insbesondere zum Ausfuhren eines 
Verfahrens nach einem der Anspriiche 1 bis 43, umfassend: 
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eine Einrichtung zum Ermitteln mittels des Verhaltensmodells, welcher 
Soll-Verlauf bestimmter Hirnaktivitatscharakteristika dem Erreichen des 
Verhaltenszieles entspricht, 

eine Einrichtung zur Erzeugung exogener elektrischer und/oder 
5 rnagnetischer Felder und Anlegen an einen Anwender mit zumindest einem 
Transmitter, 

eine Messeinrichtung zur Messung elektromagnetischer Hirnaktivitat 
mit mindestens einem Sensor, 

eine Einrichtung zur Berechnung von Hirnaktivitatscharakteristika zum 
1 0 Erreichen des Soll-Verlaufs in einem Steuer- oder Regelkreis.. 

45. Vorrichtung gemali Anspruch 44, 
dadurch gekennzeichnet, 

dass die Messeinrichtung mehrere Sensoren aufweist, die ein Sensorgitter bilden. 

15 

46. Vorrichtung nach Anspruch 44 oder 45, 
dadurch gekennzeichnet, 

dass die Einrichtung zur Erzeugung exogener elektrischer und/oder rnagnetischer 
Felder mehrere Transmitter aufweist, die ein Transmittergitter bilden. 

20 

47. Vorrichtung nach einem der Anspruche 44 bis 46, 
dadurch gekennzeichnet, 

dass mindestens eine Rechnereinheit vorgesehen ist, in der Softwareeinheiten 
zur Implementierung des Verfahrens nach einem der Anspruche 1 bis 43 
25 gespeichert sind. 

48. Vorrichtung nach einem der Anspruche 44 bis 47, 
dadurch gekennzeichnet, 

dass fur einen jeden Sensor und einen jeden Transmitter eine elektrische 
30 und/oder magnetische Abschirmung vorgesehen ist. 

49. Vorrichtung nach einem der Anspruche 44 bis 48, 
dadurch gekennzeichnet, 

dass die Messeinrichtung mechanisch von der ubrigen Vorrichtung entkoppelbar 
35 ausgebildet ist, sodass die Messeinrichtung von einem Anwender mitgefuhrt 
werden kann. 

50. Vorrichtung nach einem der Anspruche 44 bis 49, 
dadurch gekennzeichnet, 

40 dass die Sensoren und Transmitter extrakranial angeordnet sind. 
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51 . Vorrichtung nach Anspruch 50, 
dadurch gekennzeichnet, 

dass die Sensoren und Transmitter auf einer Innenseite eines die Kranialform des 
jeweiligen Anwenders nachbildenden Helmes angeordnet sind. 

5 

52. Vorrichtung nach einem der Anspruche 46 bis 51, 
dadurch gekennzeichnet, 

dass das Transmittergitter und das Sensorgitter derart verschrankt sind, dass 
jedem Transmitter Sensoren und jedem Sensor Transmitter benachbart sind. 

10 

53. Vorrichtung nach einem der Anspruche 46 bis 52, 
dadurch gekennzeichnet, 

dass im Transmittergitter Halterungen zum Aufnehmen zusatzlicher Transmitter 
vorgesehen sind, sodass die Transmitterdichte eines Transmittergitters lokal 
1 5 veranderbar ist und/oder sodass die Neigungswinkel einzelner Transmitter zum 
Kranium des Anwenders veranderbar sind. 

54. Vorrichtung nach einem der Anspruche 44 bis 53, 
dadurch gekennzeichnet, 

20 dass die Vorrichtung in eine Sensoren und Transmitter aufweisende Kopfeinheit, 
eine Zwischeneinheit und eine Basiseinheit unterteilt ist, wobei die 
Zwischeneinheit einen Rechner und Softwareeinheiten zum Ausfuhren der nicht 
kalibrierungsabhangigen Verfahrensschritte und die Basiseinheit einen Rechner 
und Softwareeinheiten zum Ausfuhren der kalibrierungsabhangigen 

25 Verfahrensschritte umfasst. 
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1. [7] AnsprQcheNr. 1-43 

well s!e slch auf Gegenstande beziehen, zu dsren Recherche die Behflrde nlcht verpflichtet 1st, namllch 

Regel 39.1(iv) PCT - Verfahren zur therapeutlschen Behandlung des 
menschlichen Oder tierischen Korpers 

2. AnsprQcheNr. 

weil sla slch auf Telle der Internatlonalen Anmeldung beziehen, die den vorgeschrlebenen Anforderungen so wenlg ent6preohen, 
daB eine sinnvolle Internationale Recherche nlcht durchgefOhrtwerden kann, namllch 



3. J^]] AnsprQcheNr. 

well es slch dabe) um abhSnglge AnsprQche handelt, die nlcht entsprechend Satz 2 und 3 der Regel 6.4 a) abgefaBt slnd. 
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" 1 — 1 zusatzllche RecherchengebOhr gerechtfertlgt hatte, hat die Behorde nlcht zur Zahlung elner solchen GebQhr aufgefordert. 



3. Da der Anmelder nur elnlge der erforderllchen zusatzilchen RecherchengebQhren rechtzeltlg entrlchtet hat, erstreckt slch 
1 — 1 Internationale Recherchenberlcht nur auf die AnsprQche, fQr die GebOhren entrlchtet worden slnd, namllch auf die 
AnsprQche Nr. 



4. [ ] Der Anmelder hat die erforderllchen zusatzilchen RecherchengebQhren nlcht rechtzeltlg entrlchtet. Der Internationale Recher- 
chenberlcht beschrankt slch daher auf die in den AnsprOchen zuerst erwahnte Erflncfcng; dlese 1st In folgenden AnsprOchen er- 
faBt: 



Bemerkungen hinslchtlich eines Wlderspruohe Q Die zusatzilchen GebOhren wurden vom Anmelder unter Wlderspruch gezahlt. 

|~] Die Zahlung zusflizllcher RecherchengebQhren erfolgte ohne Wlderspruch. 
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